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Abstract. We complete classical investigations concerning the dynamical stability of an infinite homogeneous 
gaseous medium described by the Euler-Poisson system or an infinite homogeneous stellar system described by 
the Vlasov-Poisson system (Jeans problem). To determine the stability of an infinite homogeneous stellar system 
with respect to a perturbation of wavenumber k, we apply the Nyquist method. We first consider the case of 
single-humped distributions and show that, for infinite homogeneous systems, the onset of instability is the same 
in a stellar system and in the corresponding barotropic gas, contrary to the case of inhomogeneous systems. We 
show that this result is true for any symmetric single-humped velocity distribution, not only for the Maxwellian. 
If we specialize on isothermal and polytropic distributions, analytical expressions for the growth rate, damping 
rate and pulsation period of the perturbation can be given. Then, we consider the Vlasov stability of symmetric 
and asymmetric double-humped distributions (two-stream stellar systems) and determine the stability diagrams 
depending on the degree of asymmetry. We compare these results with the Euler stability of two self-gravitating 
gaseous streams. Finally, we determine the corresponding stability diagrams in the case of plasmas and compare 
the results with self-gravitating systems. 

Key words. Stellar dynamics- hydrodynamics, instabilities 
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1. Introduction 

The dynamical stability of self-gravitating systems is an 
important problem in astrophysics. This problem was first 
considered by Jeans (1902,1929) who studied the linear 
dynamical stability of an infinite homogeneous collision- 
dominated gas described by the Euler equations. He found 
that this system becomes unstable if the wavelength of 
the perturbation exceeds a critical value A c = {-^Y^Cs 
(where c s is the velocity of sound) nowadays called the 
Jeans length. In that case, the perturbation grows ex- 
ponentially rapidly without oscillating. By contrast, for 
small wavelengths A < A c , the perturbation oscillates 
without attenuation and behaves like a gravity-modified 
sound wave. As is well-known, the Jeans approach suffers 
from a mathematical inconsistency at the start. Indeed, 
an infinite homogeneous gravitating system cannot be in 
static equilibrium since there are no pressure gradients 
to balance the gravitational force. Jeans removed this in- 
consistency by assuming that the Poisson equation de- 
scribes only the relationship between the perturbed grav- 
itational potential and the perturbed density. However, 
this assumption is ad hoc and is known as the Jeans 
swindle (Binney & Tremaine 1987). A Jeans-type anal- 



ysis can however be justified in certain situations: (i) In 
cosmology, when we work in the comoving frame, the ex- 
pansion of the universe introduces a sort of "neutralizing 
background" in the Poisson equation (like in the Jellium 
model of plasma physics). In that case, an infinite homoge- 
neous self-gravitating medium can be in static equilibrium 
(Peebles 1980). Therefore, the Jeans instability mecha- 
nism is relevant to understand the formation of galax- 
ies from the homogeneous early universe, (ii) If we con- 
sider a uniformly rotating system, the gravitational attrac- 
tion can be balanced by the centrifugal force. Therefore, 
a homogeneous system can be in static equilibrium in 
the rotating frame provided that the angular velocity and 
the density satisfy a well-defined relation (Chandrasekhar 
1961). (iii) The Jeans procedure is approximately correct 
if we only consider perturbations with small wavelengths, 
much smaller than the Jeans length (Binney & Tremaine 
1987). 

The dynamical stability of an infinite homogeneous en- 
counterless stellar system described by the Vlasov equa- 
tion was considered much later by Lynden-Bcll (1962)0 
who used methods similar to those introduced by Landau 



See also Simon (1961) and Liboff (1963). 
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(1946) in plasma physics. There exists indeed many analo- 
gies between self-gravitating systems and plasmas since 
the Coulombian and Newtonian interactions both corre- 
spond to an inverse square law. However, there exists 
also crucial differences. First of all, gravity is attractive 
whereas electricity is repulsive. On the other hand, as al- 
ready indicated, the self-gravity of a uniform gravitating 
system is usually not neutralized, contrary to a plasma 
where the presence of electrons and ions provide electri- 
cal neutrality. In order to circumvent the Jeans swindle, 
Lynden-Bell considered the case of a rotating system and 
noted that the angular velocity plays a role similar to 
the magnetic field in plasma physics (Bernstein 1958). He 
found that the system becomes unstable above a criti- 
cal wavelength. For a Maxwellian distribution, this criti- 
cal wavelength A c = (^) 1 ^ 2 c is equivalent to the Jeans 
length if we identify the r.m.s. velocity of the stars a in one 
direction with the velocity of sound c s in an isothermal 
gas. In that case, the perturbation grows exponentially 
rapidly without oscillating. By contrast, for small wave- 
lengths A < A c , the perturbation is damped exponentially. 
This is similar to the Landau damping in plasma physics 
that is derived in the complete absence of collisions. 

Sweet (1963) considered the gravitational instability of 
a system of gas and stars in relative motion and provided 
a detailed analysis of the dispersion relation in different 
cases. He considered in particular the situation where the 
gas is at rest (U — 0) and the star system is made of 
two identical interpenetrating streams with Maxwellian 
distribution moving in opposite direction with equal ve- 
locity ±v a . This corresponds to the Kapteyn-Eddington 
two-stream representation in the solar neighborhood. He 
mentioned that the two humps could be asymmetric but 
he did not study the consequence of this asymmetry in de- 
tail. One important conclusion of his work is that the crit- 
ical Jeans wavelength is reduced by the relative motion of 
the gas and stars and that it vanishes when the relative ve- 
locity v a exceeds a limit of the order of the effective veloc- 
ity of sound c s in the gas. In particular, the star streaming 
in the solar neighborhood can cause the gas to be unsta- 
ble at all wavelengths. A similar conclusion was reached by 
Talwar & Kalra (1966) who considered the contrastream- 
ing instability of two self-gravitating gaseous streams with 
velocity ±U described by fluid equations. In the subsonic 
regime U < c s , they found that the critical Jeans wave- 
length is reduced and tends to zero as the streaming ve- 
locity U approaches the velocity of sound c s (by contrast, 
in the supersonic regime U > c s , the critical wavelength 
is larger than the Jeans length without streaming). The 
works of Sweet (1963) and Talwar & Kalra (1966) were 
further developed by Ikeuchi ct al. (1974) who presented 
various stability diagrams for two-stream stellar systems, 
gaseous-stellar systems, two gaseous systems and plasmas. 
For a purely two-stream stellar systems (without gas), 
they found that the critical Jeans length becomes larger 
due to the stabilization effect of relative velocity, contrary 
to the case where a gas component is present. The two- 
stream instability was also discussed by Araki (1987) for 



infinite homogeneous and uniformly rotating stellar sys- 
tems. 

The seminal works of Jeans (for gaseous systems) and 
Lynden-Bell (for stellar systems) have been continued 
in several directions. For example, the Jeans instability 
of a self-gravitating gas in the presence of a magnetic 
field or an overall rotation has been studied in detail by 
Chandrasckhar (1961). On the other hand, the gravita- 
tional stability of a disk of stars is treated in the classical 
paper of Toomre (1964). Other interesting contributions 
are reviewed by Binncy & Trcmainc (1987). However, this 
problem is largely academic because real stellar systems 
(like stars and galaxies) are not spatially homogeneous 
and are limited in space. Now, when spatial inhomogene- 
ity is properly taken into account, the picture becomes 
very different. The dynamical stability of stars has been 
considered by various authors such as Eddington (1918), 
Ledoux (1945) and Chandrasekhar (1964) and different 
stability criteria for the Euler-Poisson system have been 
obtained. For example, it can be shown that polytropic 
stars are Euler stable if their index n < 3 while they arc 
unstable if 3 < n < 5 (polytropic stars with n > 5, in- 
cluding isothermal stars n — > +oo, have infinite mass). 
On the other hand, the Vlasov stability of stellar sys- 
tems has been studied by Antonov (1960a), Doremus et al. 
(1971), Kandrup & Sygnet (1985) and others. They have 
shown that all stellar systems with a distribution function 
/ = /(e) which is a monotonically decreasing function of 
the energy are linearly dynamically stable with respect 
to the Vlasov-Poisson system. In particular, all polytropic 
galaxies with index in the range 3/2 < n < 5 (such that 
the total mass is finite) arc stablqj. Very recently, it was 
proven rigorously by Lemou et al. (2009) that these sys- 
tems are also nonlinearly dynamically stable with respect 
to the Vlasov-Poisson system. These results show that, 
when the spatial inhomogeneity of the system is properly 
accounted for, the Jeans instability is suppressed^. 

Recently, there was a renewed interest for this classi- 
cal problem due to analogies with other physical systems. 
For example, the community of statistical mechanics in- 
volved in the dynamics and thermodynamics of systems 
with long-range interactions (Dauxois et al. 2002) has 
studied a toy model called the Hamiltonian Mean Field 
(HMF) model which presents many features similar to self- 
gravitating system^. In particular, this model presents an 
instability below a critical temperature that is similar to 

2 The difference between the dynamical stability of gaseous 
stars and stellar systems, which is related to the Antonov 
(1960b) first law, can also be related to a notion of "ensemble 
inequivalence" like in thermodynamics (see Chavanis 2006a). 

' 5 A form of Jeans instability is recovered for spatially inho- 
mogeneous isothermal (n = oo) and polytropic (n > 5) systems 
confined within a box (Padmanabhan f 990, Semelin et al. 200f , 
Chavanis 2002a,2002b, Taruya & Sakagami 2003a). 

4 In fact, this model is directly inspired by astrophysics. It 
was introduced by Pichon (1994) as a simplified model to de- 
scribe the formation of bars in disk galaxies (see a detailed 
historic in Chavanis et al. (2005)). 
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the Jeans instability in astrophysics. Interestingly, a spa- 
tially homogeneous phase exists for this model at any tem- 
perature (it is stable for T > T c and unstable for T < T c ). 
Therefore, there is no mathematical inconsistency (i.e. no 
"Jeans swindle" ) when we perform the stability analysis of 
the homogeneous phase of this system (Inagaki & Konishi 
1993, Choi & Choi 2003, Yamaguchi ct al. 2004, Chavanis 
& Delfini 2009). On the other hand, in biology, several re- 
searchers in physics and applied mathematics have stud- 
ied the phenomenon of chemotaxis using the Kcller-Segcl 
(1970) model in order to explain the self-organization of 
bacterial colonies, cells or even social insects. Recently, 
it was shown that the chemotactic instability in biology 
bears some deep analogies with the Jeans instability in 
astrophysics (Chavanis 2006b Interestingly, in the bio- 
logical problem, there is no "Jeans swindle" because the 
Poisson equation is replaced by a more complex reaction- 
diffusion equation that allows for the existence of infinite 
and homogeneous distributions. 

These analogies were a source of motivation to study 
anew the classical Jeans problem in astrophysics (for the 
Eulcr-Poisson and Vlasov- Poisson systems). In fact, we 
discovered that the study of this classical problem was 
still incomplete. For example, the case of (spatially homo- 
geneous) polytropic distribution functions has not been 
investigated in detail and the case of a symmetric double- 
humped distribution considered by Ikeuchi et al. (1974) 
needs further discussion. On the other hand, the stability 
of an asymmetric double-humped distribution has never 
been considered in the astrophysical literature (although 
the interest of such distributions was mentioned early by 
Sweet (1963)) and could form an interesting theoretical 
problem. The stability criteria for homogeneous stellar 
systems can be established very easily with the powerful 
Nyquist method (Nyquist 1932) used in plasma physics. 
This method can be readily applied to homogeneous self- 
gravitating systems. However, since the gravitational in- 
teraction is attractive while the electric interaction is re- 
pulsive, a crucial sign changes in the dispersion relation 
and the problem must be reconsidered. In particular, this 
change of sign is the reason for the Jeans instability in as- 
trophysics and the Nyquist method gives a nice graphical 
illustration of this instability. This method does not seem 
to be well-known by astrophysicists and is not exposed in 
classical monograph^. Although it essentially has an aca- 
demic interest since real stellar systems are not spatially 
homogeneous^, we think that it deserves a detailed peda- 

5 Curiously, Keller & Segel (1970) did not point out this 
analogy and interpreted instead the chemotactic instability in 
relation to Benard convection. 

6 As mentioned by the referee, an early application of the 
Nyquist method was made by Lynden-Bell (1967) in a not very 
easily accessible paper. Nyquist's method has also been used 
by Weinberg (1991) to study the stability of spherical stellar 
systems. We are not aware of other references. 

7 In fact, the Nyquist method can be generalized to apply to 
spatially inhomogeneous self-gravitating systems (see Weinberg 
1991). 



gogical exposition. We thus propose an exhaustive study 
of the linear dynamical stability of infinite and homoge- 
neous gaseous (Euler) and stellar (Vlasov) systems that 
uses the Nyquist method and completes previous works 
on that problem. 

The paper is organized as follows. In Sec. [5J we con- 
sider a self-gravitating barotropic gas described by the 
Eulcr-Poisson system. We recall the classical Jeans insta- 
bility criterion and, for future comparison, consider ex- 
plicitly the case of isothermal and polytropic equations 
of state. In Sec. 131 we consider a collisionless stellar sys- 
tem described by the Vlasov-Poisson system. Using the 
Nyquist method, we derive the Jeans instability criterion 
for an infinite and homogeneous distribution. We show 
that, for any single humped distribution function of the 
form / = f(v 2 ) with f'{v 2 ) < 0, the Jeans instability cri- 
terion for a stellar system is equivalent to the Jeans insta- 
bility criterion for the corresponding barotropic gas with 
the same equation of state. Of course, the dispersion rela- 
tion and the evolution of the perturbation is different in 
the two systems (gas and stellar system) but the threshold 
of the instability is the same. This generalizes the result 
obtained by Lynden-Bell (1962) for the Maxwellian dis- 
tribution. In Sees. 2] and [SJ we specifically consider the 
case of isothermal and polytropic distribution functions 
and derive analytical expressions for the growth rate and 
damping rate of the perturbation by adapting methods of 
plasma physics to the gravitational context. In Sec. [5J we 
consider a symmetric double-humped distribution made 
of two counter-streaming Maxwellians with velocities ±v a 
and establish the stability diagram. In Sec. we generalize 
our results to the case of an asymmetric double-humped 
distribution. We find that: (i) The system is stable for 
A < (A c )* and unstable for A > (A c )* where the critical 
wavelength (Ac)* is a non-trivial function of the relative 
velocity v a and asymmetry parameter A. (ii) The critical 
wavelength (A c )* is always larger than the critical wave- 
length in the absence of streaming (v a = 0) so that the in- 
stability is delayed, (iii) The phase velocity of the unstable 
mode corresponds to the global maximum of the distribu- 
tion function. In Sec. [51 we consider the case of plasmas. 
We recall the well-known fact that single-humped distri- 
butions are always stable. Then, we consider the case of 
symmetric and asymmetric double-humped distributions. 
For an asymmetry parameter A < A* = 3.3105: (i) the 
system is stable for v 2 < y c (A)T where y c (A) depends on 
the asymmetry (y c = 1.708 for a symmetric distribution); 
(ii) for v 2 > y c (A)T, the system is stable for A < (A c )* 
and unstable for A > (A c )* where (A c )* depends on the 
relative velocity v a and on the asymmetry parameter A. 
For an asymmetry parameter A > A* = 3.3105: (i) the 
system is stable for v 2 < y*(A)T where y*(A) depends on 
the asymmetry; (ii) for y*(A)T < v 2 < y c (A)T, the sys- 
tem is stable for A < (A c )i, unstable for (A c )i < A < (A c )2 
and stable for A > (A c )2 where (A c )i and (A c )2 depend on 
the relative velocity v a and on the asymmetry parameter 
A; (iii) for v 2 > y c (A)T, the system is stable for A < (A c )i 
and unstable for A < (A c )i- Throughout the paper, we also 
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compare our results obtained for stellar systems (Vlasov) 
to those obtained for gaseous systems (Euler). 

2. Self-gravitating barotropic gas 

2.1. The Euler-Poisson system 

We consider a self-gravitating gas described by the Euler- 
Poisson system 



+ V • {pa) = 0, 



dp 
&t 

du , . 1 

— + u-V u = — Vp-V$, 

at p 

A$ = AirGp. 



(1) 

(2) 
(3) 



These are the usual equations used to study the dynamical 
stability of a gaseous medium like a molecular cloud or 
a star. We assume that the gas is a perfect gas in local 
thermodynamic equilibrium (L.T.E). Therefore, at each 
point, the distribution function of the particles (atoms or 
molecules) is of the form 



/(r,v,i) 



1 



2-KT{v,t) 



3/2 



p(r,t)e~ 



^— »(r,t)]' ! 
2T(r,t) 



(4) 



Note that the kinetic temperature T(yA) = -|(w 2 ) where 
w = v — u(r, t) is the relative velocity. Introducing the 
pressure P — \ j fw 2 dv, we find that the local equation 



of state is 
p(r,t) = p(r,t)T(r,t). 



(5) 



The system of equations ([I])-© is not closed since it must 
be supplemented by an equation for the energy or, equiv- 
alently, for the temperature T(r, t). Here, we shall restrict 
ourselves to the case of a barotropic gas for which the pres- 
sure depends only on the density 



=p[p( r ,t)]. 

The local velocity of sound is 
c 2 (r,t) =p'(p(r,t)). 

A steady state of the Euler-Poisson system satisfies u 
and the condition of hydrostatic equilibrium 

Vp + pV$ = 0. 



(6) 



(7) 



(8) 



Combining the condition of hydrostatic equilibrium ((5]) 
and the equation of state p = p(p), we find that the density 



8 Throughout this paper, by an abuse of language, we shall 
define the kinetic temperature by T(r,t) = 4((v — u(r, t)) 2 } 
where d is the dimension of space. Therefore, the kinetic tem- 
perature is a measure of the local velocity dispersion of the 
particles in one dimension. It is related to the real kinetic tem- 
perature by T(r,t) — kBT rea i(r,t)/m. We do that in order to 
extend this definition to the case of collisionless stellar systems 
described by the Vlasov equation where the individual mass of 
the stars does not appear. 



is a function p = p(&) of the gravitational potential given 

by 
rPp'(p') 



P' 



-dp' 



(9) 



We note that p'{p) = -p/p'($). Since p'{p) = c 2 > 0, we 
find that //($) < 0. The density is a monotonically de- 
creasing function of the gravitational potential. The total 
energy of the star is 

W[p,n] = J ' p J* ?^ldp'dv+ l - J p^dr + J p^-dr, 

(10) 

including the internal energy, the potential energy and 
the kinetic energy of the mean motion. The mass M[p] 
and the energy W[p, u] are conserved by the barotropic 
Euler-Poisson system: M — W — 0. Therefore, a mini- 
mum of the energy functional at fixed mass determines a 
steady state of the barotropic Euler-Poisson system that 
is nonlinearly dynamically stable. Here, we have stability 
in the sense of Lyapunov. This means that the size of the 
perturbation is bounded by the size of the initial pertur- 
bation for all times. We are led therefore to considering 
the minimization problem 



min 



M[p] = M}. 



(11) 



The critical points of energy at fixed mass are determined 
by the Euler-Lagrange equation <5W — pSM = where p 
is a Lagrange multiplier. This yields the condition of hy- 
drostatic equilibrium . Then, a critical point of energy 
at fixed mass is a (local) energy minimum iff 

5 2 W= J ^-(5p) 2 dr+^ J 6p6$dr>0 (12) 

for all perturbations Sp that conserve mass: SM = 0. 

2.2. The Jeans instability criterion 

We consider the linear dynamical stability of an infinite 
homogeneous gaseous medium described by the hydrody- 
namic equations ([1]), ([2]), ((3]) and . This is the clas- 
sical Jeans (1902,1929) problem. Linearizing Eqs. (P)-© 
around a solution p = $ = const, and u = 0, making the 
Jeans swindle and decomposing the perturbations in nor- 
mal modes of the form e l ( k r - w t) ; we obtain the classical 
dispersion relation (Binney & Tremaine 1987): 

u 2 = c 2 s k 2 - AuGp, (13) 

where c 2 = dp/ dp is the velocity of sound. Since lu 2 is real, 
the complex pulsation tu is either real or purely imaginary. 
The condition uj = determines a critical wavenumber 

kl - (14) 

called the Jeans wavenumber for a barotropic gas. The 
dispersion relation can be rewritten 

uj 2 _ k 2 _ 

4rrGp " fc c 2 ~ (15) 
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For k > k c , we have to = co r = ±\/ljj 2 so that the pertur- 
bation oscillates like e ~ lu)rt with a pulsation 



±^c 2 s k 2 - AnGp, 



(16) 



without attenuation. This corresponds to gravity-modified 
sound waves. In that case the system is stable. For k < k c , 



wc have to = l 
like e Uit with a growth rate 

LOi = y 4nGp — c 2 k 2 , 



so that the perturbation grows 



(17) 



without oscillation (the second mode is attenuated expo- 
nentially rapidly so that the growing mode dominates). 
In that case, the system is unstable. This is the so-called 
Jeans instability. The growth rate is maximum for k = 
and its value is Wi(k = 0) = y/AitGp. 

In conclusion, a perturbation with wavenumber k is 
stable if 



(18) 



and unstable otherwise. These results are valid for an ar- 
bitrary barotropic equation of state p = p{p). We now 
specialize on particular cases, namely isothermal and poly- 
tropic equations of state. 

2.3. Isothermal equation of state 

We first consider an isothermal equation of state 

p(r,t) =p(r,t)T, (19) 

where the temperature T(r, t) — T is a uniform. At hydro- 
static equilibrium, according to Eq. ([5]), the relation be- 
tween the density and the gravitational potential is given 
by the Boltzmann law 

p = A'e-*, (20) 
where A' is a constant. The energy functional (|10p reads 

W[p,u] = k B T J plnprfr+i / p<S>dr + / p\ dv , 



and it can be viewed as a Boltzmann free energy F = 
E — TS where E is the macroscopic energy and S the 
Boltzmann entropy. The velocity of sound (JJJ is uniform: 



T. 



(22) 



Finally, for an infinite homogeneous isothermal gas, the 
Jeans instability criterion takes the classical form 



e < k\ = ^e. 

C rjn 



(23) 



2.4. Polytropic equation of state 

The equation of state of a polytropic gas@ is 
p{v,t)=Kp{v,t)\ 



(24) 



where K is a constant. The polytropic index n is defined 
by 7 = 1 + 1/n. For n — > +oo, we recover the isothermal 
case with 7 = 1. At hydrostatic equilibrium, according to 
Eq. ©, the relation between the density and the gravita- 
tional potential can be written 



A 



7 



1 



(25) 



where A is a constant. The energy functional (|10j) can be 
put in the form 



W[p,u] 



K 
~ 1 



{p 1 - p)dr + 



p§dr 



p—dr. 



(26) 



In the limit 7 — > 1, we recover the results (j2"0"]) and (|2"Tj) 
obtained for isothermal distributions. For a polytropic gas, 
the local temperature T(r, t) given by Eq. © reads 



T(v,t)=Kp{v,t) 



7-1 



(27) 



We note that the temperature is position dependent (while 
the specific entropy s is uniform) and related to the den- 
sity by a power law (this is the local version of the usual 
isentropic law TV~ l = const, in thermodynamics). The 
polytropic index n is related to the gradients of tempera- 
ture and density according to 



VlnT = -Vlnp. 

n 



(28) 



Combining Eqs. (|2"5j) and ([2"7]l. we note that, for a poly- 
tropic gas at equilibrium, the relation between the kinetic 
temperature and the gravitational potential is linear so 
that 



7 - 1 
VT= V$. 

7 



(29) 



The coefficient of proportionality is related to the adia- 
batic index 7. The velocity of sound ([7]) can be expressed 



(21) c s (r,t) 2 =jT(r,t), 



(30) 



9 The polytropic equation of state corresponds to an adia- 
batic transformation for which the specific entropy s(r,t) = s 
is uniform. In that case, 7 = c p /c v is the ratio of specific heats 
at constant pressure and constant volume. For a monoatomic 
gas 7 = 5/3. This law describes convective regions of stellar in- 
terior. An approximate polytropic equation of state with index 
7 ~ 3.25 also holds in the radiative region of a star. Finally, 
classical white dwarf stars are described by a polytropic equa- 
tion of state with index 7 = 5/3 (n = 3/2). This equation of 
state is valid for a completely degenerate gas of electrons de- 
scribed by the Fermi-Dirac statistics at T = (Chandrasekhar 
1942). 
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and it is usually position dependent. However, when we 
consider an infinite and homogeneous distribution (mak- 
ing the Jeans swindle) , the velocity of sound and the tem- 
perature are uniform. In that case, we can speak of "the 
temperature T of the polytropic gas" and write the Jeans 
instability criterion as 



k z <ki = 



7T 



(31) 



We thus find that the critical Jeans wavenumber kc of 
a polytropic gas is related to the critical Jeans wavenum- 
ber kc ls °^ of an isothermal gas with the same temperature 
Thy 

k (poly) = _fc(<»°) ) (32) 

where the proportionality factor is the polytropic index 7. 
The positivity of c 2 = dp/ dp implies that 7 > 0, i.e. n > 
or n < —1 (for — 1 < n < 0, the gas is always unstable, 
even in the absence of gravity, since ui 2 < 0). For 7 = 
(n = —1), the gas is unstable to all wavelengths since 
ki polv) = +00. For < 7 < 1 (n < -1), k { c poly) > k { c zso) . 
For 7 = 1 (n = 00), kc P ° ly ' > = k^ ^ 1 (isothermal case). For 
7 > 1 (n > 0), ki poly) < ki zso) . For 7 = +00 (n = 0), 
the gas is stable to all wavelengths since fcc P ° iy - ) = (solid 
medium). For an isentropic gas for which 7 = c p /c v , the 
Mayer relation c p — c v = ks implies that 7 > 1 (n > 0). 
Therefore, for an isentropic gas, the gravitational instabil- 
ity is always delayed (i.e., it occurs for larger wavelengths) 
with respect to an isothermal gas with the same temper- 
ature. 

3. Collisionless stellar systems 

3.1. The Vlasov-Poisson system 

We consider a self-gravitating system described by the 
Vlasov-Poisson system 

df df df 
— + v • — + F • — =0, 
dt dr dv 



A$ = 4ttG / fdv 



(33) 



(34) 



where F = — V<I> is the self-consistent gravitational field 
produced by the particles. The Vlasov description assumes 
that the evolution of the system is encounterless. This is a 
very good approximation for most astrophysical systems 
like galaxies and dark matter because the relaxation time 
due to close encounters is in general larger than the age 
of the universe by several orders of magnitude (Binney & 
Trcmainc 1987). The Vlasov equation admits an infinite 
number of stationary solutions given by the Jeans theo- 
rem (Jeans 1915). For example, distribution functions of 
the form / = /(e) which depend only on the individual 
energy e = v 2 /2 + 3>(r) of the stars are particular steady 
states of the Vlasov equation. They describe spherical stel- 
lar systems. 



The Vlasov equation conserves the mass M[f] = 
J f drdv, the energy E[f] = | / fv 2 drdv + i J p<& dr and 
the Casimirs //, = J h(f) drdv for any continuous function 
h. Let us introduce the functionals 



C(f) drdv, 



(35) 



where C is a convex function (C" > 0). These particular 
Casimirs will be called "pseudo entropies'^- 

The two constraints problem: since the Vlasov equation 
conserves M, E and S, the maximization problem 

max{5[/] I E[f] = E, M[f] = M} , (36) 

determines a steady state of the Vlasov-Poisson system 
that is nonlinearly dynamically stable. The critical points 
of pseudo entropy at fixed mass and energy are deter- 
mined by the variational principle SS — /35E — aSM = 0, 
where (3 (pseudo inverse temperature) and a (pseudo 
chemical potential) are Lagrange multipliers. This yields 
C'(f) = —f3e — a. Since C is convex, this relation can be 
reversed to give / = F(f3e+a) where F = (C') _1 (— x). We 
have /'(e) = -f}/C"(f). Therefore, if P > 0, the critical 
points of pseudo entropy at fixed mass and energy deter- 
mine distributions of the form / = /(e) with /'(e) < 0: the 
distribution function is a monotonically decreasing func- 
tion of the energjO Then, a critical point of pseudo en- 
tropy at fixed mass and energy is a (local) maximum iff 

52 S = -\ J C"(f)(Sf) 2 drdv-±p [s P 6$dr<0, (37) 

for all perturbations Sf that conserve mass and energy at 
first order: SM = SE = 0. 

The one constraint problem: the minimization problem 



min {-F[/] 



E[f]-TS[f] I M[f]=M}. 



(38) 



where T = 1/p > is prescribed, determines a steady 
state of the Vlasov-Poisson system that is nonlinearly dy- 
namically stable. The functional F[f] is a pseudo free en- 
ergy. The critical points of pseudo free energy at fixed mass 
are determined by the variational principle 5F + aTSM = 
0, where a (pseudo chemical potential) is a Lagrange mul- 
tiplier. This yields C'(f) = —Pe — a. Then, a critical point 
of pseudo free energy at fixed mass is a (local) minimum 
iff 

S 2 F =±tJ C"(f)(5f) 2 drdv + l - J 5pSf> dr > 0, (39) 
for all perturbations 5f that conserve mass: SM = 0. 



10 The functionals H = —JC(f)drdv defined in terms of 
the coarse-grained distribution / are called generalized H- 
functions (Tremaine et al. 1986). 

11 More generally, the solutions of (|36p are of the form / = 
/(e) where / is monotonic, decreasing at positive temperatures 
and increasing at negative temperatures. For realistic stellar 
systems, the DF should decrease close to the escape energy 
e = 0. Therefore, for the class of distributions considered, / 
must decrease everywhere implying ft > 0. 
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The no constraint problem: the maximization problem 



max {G[f} = S[f] - pE[f\ ~ ^M[f]} , 



(40) 



where j3 and a are prescribed, determines a steady state of 
the Vlasov-Poisson system that is nonlinearly dynamically 
stable. The functional G[f] is a pseudo grand potential. 
The critical points of pseudo grand potential G, satisfying 
8G = 0, are given by C'(f) = —fie — a. Then, a critical 
point of pseudo grand potential is a (local) maximum iff 



5 2 G = -- 



\ J C"(f)(5f) 2 drdv-^f3 J 5p5$dv<Q, (41) 



for all perturbations 5f. 

The optimization problems (|38p and (|40|l have the 
same critical points (canceling the first order variations). 
Furthermore, a solution of (j4"0"]) is always a solution of the 
more constrained dual problem (|38[) . Indeed, if inequality 
pTT) is true for all perturbations, it is true a fortiori for 
all perturbations that conserve mass. Similarly, a solution 
of (|3"5)) is always a solution of the more constrained dual 
problem (|36[) . Indeed, if inequality (|39[) is true for all per- 
turbations that conserve mass, it is true a fortiori for all 
perturbations that conserve mass and energy. However, 
the reciprocal is wrong. A solution of (|3"5)) is not neces- 
sarily a solution of (|4T)]) , and a solution of (|36|) is not nec- 
essarily a solution of (|38[), This is similar to the notion 
of ensemble inequivalence in thermodynamics (Ellis et al. 
2000, Bouchct & Barre 2005, Chavanis 2006a). Indeed, the 
two constraints problem (|36[) is similar to a condition of 
microcanonical stability, the one constraint problem ()38l) 
is similar to a condition of canonical stability, and the no 
constraint problem (|40[) is similar to a condition of grand 
canonical stability. The implication (gOJ) => => (EH) 
is similar to the fact that grand canonical stability im- 
plies canonical stability which itself implies microcanon- 
ical stability (but not the converse) in thermodynamics. 
Therefore, (|40[) provides just a sufficient condition of non- 
linear dynamical stability that is less refined than (|38[) . 
and (|38[) provides just a sufficient condition of nonlinear 
dynamical stability that is less refined than (1361) . 

The most refined problem: the minimization problem 



mm 

/ 



{E[f] | all the Casimirs Ih} , 



(42) 



determines a steady state of the Vlasov-Poisson sys- 
tem that is nonlinearly dynamically stable. A distribu- 
tion function is a critical point of energy for symplec- 
tic perturbations (i.e. perturbations that conserve all the 
Casimirs) iff /(r,v) is a steady state of the Vlasov equa- 
tion (Bartholomew 1971, Kandrup 1991). Furthermore, if 
we consider spherical stellar systems for which / = /(e), 
it can be shown (Bartholomew 1971, Kandrup 1991) that 
a critical point of energy for symplectic perturbations is a 
(local) minimum iff 



5 Z E = - 



(Sf) 2 



dvdw 



5p5$dr > 0, 



(43) 



for all perturbations 5f that conserve energy and all the 
Casimirs at first order: 5E = SIh = 0. Such symplec- 
tic (physically accessible) perturbations are of the form 
Sf = D5g where D = v • V r — V<I> • V v is the advective 
operator in phase space and Sg(r, v) is any perturbation. 
Inequality (|4"3")l corresponds to the Antonov (1960) crite- 
rion of dynamical stability that was obtained by investi- 
gating the linear dynamical stability of a steady state of 
the Vlasov equation. Since this is the most constrained cri- 
terion, this is the most refined one. We note that inequal- 
ity (f4*3")l is equivalent to inequalities (|3"T|) . (|3"9"]l and f4"Tj) 
if we use the identity /'(e) = —/3/C"(f) derived above. 
However, the classes of perturbations to consider are dif- 
ferent: in (jUJ), we must consider all perturbations, in (|39l 
we must consider perturbations that conserve mass, in 
(|37[) we must consider perturbations that conserve mass 
and energy, and in (|43[) we must consider perturbations 
that conserve energy and all the Casimirs. Therefore, we 
have the implications (|ID )) =»• ([55 ]) =4> (|56 ) =»• (|32" |) . The 
connection between (|4"2"|) and the preceding optimization 
problems can be understood as follows. The minimization 
problem (|4"2"|) . conserving all the Casimirs, is clearly more 
refined than the minimization problem 



min {£[/] 



M[f]=M, S[f}=S}. 



(44) 



where only the mass and one Casimir of the form Q35[) is 
conserved. Now, it is easy to show that the minimization 
problem (|44[) is equivalent to the maximization problem 
(EU). Hence, we have the chain of relations 



lHB=>!|3HH=>®^!@m=* 422)- 



(45) 



To summarize, the minimization problem (|42p is the most 
refined stability criterion because it tells that, in order to 
settle the dynamical stability of a stellar system, we just 
need considering symplectic (i.e. dynamically accessible) 
perturbations. Of course, if inequality (|43[) is satisfied by a 
larger class of perturbations, as implied by problems (I40p , 
(|38p. (|3l)|) and (|4~4"]) . the system will be stable a fortiori. 
Therefore, we have the implications (|4"5]) . Problems (|4"0)) . 
(|3"5)h (|3"6")) and (|4"4")l provide sufficient (but not necessary) 
conditions of dynamical stability. A steady state can be 
stable according to (f3"5)) while it does not satisfy (|4T)]) , or 
it can be stable according to ([3"6"]l and (|44p while it does not 
satisfy (|3"5)) . or it can be stable according to (|4"2"j) while it 
does not satisfy (f3l>|) and ([4"4")l . Therefore, the criterion (|42p 
is more refined than ([44]) or ((36)) . which is itself more re- 
fined than (f38)) . which is itself more refined than (j40]) . Let 
us give an astrophysical illustration of all that (Chavanis 
2006a). Using (|3"5)) . we can show that stellar polytropes 
with 3/2 < n < 3 arc dynamically stable because they are 
minima of F at fixed mass M. By contrast, stellar poly- 
tropes with 3 < n < 5 are not minima of F at fixed mass. 
However, using (|3"6")) . we can show that stellar polytropes 
with 3/2 < n < 5 are dynamically stable because they 
are maxima of S at fixed mass and energy. However, not 
all stellar systems of the form / = /(e) are maxima of S 
at fixed mass and energy. But, using (|42|). we can show 
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that all spherical galaxies / = /(e) with /'(e) < are dy- 
namically stable. This last statement has been shown for 
linear stability by Kandrup (1991). A rigorous mathemat- 
ical proof has been given recently by Lemou et al. (2009) 
for nonlinear stability. 

Remark: similar results are obtained in 2D fluid me- 
chanics based on the Euler-Poisson system (see Chavanis 
2009). Criterion (j4"2"]l is equivalent to the so-called Kelvin- 
Arnol'd energy principle, criterion (|40[) is equivalent to the 
standard Casimir-energy method (see Holm et al. 1985) 
introduced by Arnol'd (1966) and criterion ((55]) is equiv- 
alent to the refined stability criterion given by Ellis et al. 
(2002). 

3.2. The corresponding barotropic star 

To any stellar system with / = /(e) and /'(e) < 0, we can 
associate a corresponding barotropic star with the same 
equilibrium density distribution (Lynden-Bell & Sanitt 
1969). Indeed, defining the density and the pressure by 
p = J fdw and p = g J fv 2 dv, we get p = p($) and 
p = p($). Eliminating the potential $ between these two 
expressions, we find that p — p{p). Furthermore, taking 
the gradient of the pressure and using the chain of identi- 
ties 



-V$ 



9/ 
dv 



v e?v 



f'(e)v 2 dv 



(46) 



we obtain the condition of hydrostatic equilibrium ([8]). 
Finally, we define the kinetic temperature of an isotropic 
stellar system by the relation 



Therefore, the quantity 



T(r>4 



Ufv^dv p(r) 



p(r)' 



(47) 



(48) 



measures the velocity dispersion (in one direction) of an 
isotropic stellar system. In general, the kinetic tempera- 
ture is position dependent. 

Finally, we can show that the variational principles 
(|38l) and (fTTj) are equivalent, i.e. a stellar system is a min- 
imum of pseudo free energy at fixed mass iff the corre- 
sponding barotropic gas is a minimum of energy at fixed 
mass (see Appendix . This leads to a nonlinear version 
of the Antonov (1960b) first law: "a stellar system with 
/ = /(e) and /'(e) < is nonlinearly dynamically stable 
with respect to the Vlasov-Poisson system if the corre- 
sponding barotropic gas is nonlinearly dynamically stable 
with respect to the Euler-Poisson system" . However, the 
reciprocal is wrong because, as we have already indicated, 
(|38[) provides just a sufficient condition of nonlinear dy- 
namical stability with respect to the Vlasov equation. A 
galaxy can be dynamically stable according to criterion 



(l36l) [or even more generally (|4"2"j) ] while it fails to satisfy 
criterion ([55]) . Therefore, the nonlinear Antonov first law 
is similar to a notion of ensembles inequivalence between 
microcanonical and canonical ensembles in thermodynam- 
ics (Chavanis 2006a). 

3.3. The dispersion relation 

We shall study the linear dynamical stability of a spatially 
homogeneous stationary solution of the Vlasov equation 
described by a distribution function / = /(v). Like for 
an infinite homogeneous gas, we make the Jeans swindle. 
Linearizing the Vlasov equation around this steady state, 
taking the Laplace-Fourier transform of this equation and 
writing the perturbations in the form gHk- 1 "-"*)^ we ob- 
tain the classical dispersion relation (Binney & Trcmainc 
1987): 



e(k,w) = 1 



AttG 



k- ^ 



-dv = 0, 



(49) 



where e(k,w) is similar to the "dielectric function" of 
plasma physics (with the sign + instead of — ). For a given 
distribution /(v), this equation determines the complex 
pulsation(s) uj = ui r + iuii of a perturbation with wavevec- 
tor k. Since the time dependence of the perturbation is 
Sf ~ e~ tUrt e^ it , the system is linearly stable if uii < 
and linearly unstable if u>i > 0. The condition of marginal 
stability corresponds to uii = 0. 

If we take the wavevector k along the z-axif0, the 
dispersion relation becomes 



AttG f k§f 

e(k,Lo) = 1 + —r- I ^^dv x dv y dv z = 0, 

k J c kv z — uj 



(50) 



where the integral must be performed along the Landau 
contour C (Binney & Tremaine 1987). We define 



g(v z ) = / fdv x dv y . 



(51) 



In the following, we shall note v instead of v z and / instead 
of g. With these conventions, the dispersion relation (|5Q|) 
can be rewritten 



e{k,uj) = 1 + -p- / —dv = 



G 



(52) 



This is the fundamental equation of the problem. For fu- 
ture reference, let us recall that 



dv 



/» 



dv, {u>i > 0), 



(53) 



f'(v) r+ °° *' 

±^rdv = P 



/'(«) 



-00 " k 



dv + iirf , {u>i = 0), 



(54) 



12 If the distribution function is isotropic, there is no restric- 
tion in making this choice. 
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<) 



dv 



cfo + 2ttz/' 



(w< < 0), 



(55) 

where P denotes the principal value. 

In general, the dispersion relation (f5"2"j) cannot be 
solved explicitly to obtain uj(k) except in some very sim- 
ple cases [3 For example, for cold systems described by 
the distribution function / = pS(v — Vo), we obtain after 
an integration by parts 

to = vok ± iy/ AirGp. (56) 
In particular, when vq = 0, we get 



to = ±i^AnGp. (57) 

The system is unstable to all wavenumbcrs and the per- 
turbation grows with a growth rate u>i = \JknGp. We 
also note that the dispersion relation ([57)1 coincides with 
the dispersion relation (fT7| of a self-gravitating gas with 
c s = 0. 

3.4. The condition of marginal stability 

For uji = 0, the real and imaginary parts of the dielectric 
function e(k,u> r ) = e r (k,uj r ) + iei{k,u> r ) are given by 



4ttG 

e r (k,uj r ) = 1 + —r^-P 



ei(k,oj r ) = ^ f (u) r /k) 



m 

v — oj r /k 



dv, 



(58) 



(59) 



The condition of marginal stability corresponds to 
e(k, bS) = and = 0, i.e. e r (k,uj r ) — Ei(k,uj r ) = 0. 
We obtain therefore the equations 



4ttG f + °° f'(v) , 
l + —^P I - V ; ,. dv = 0, 



k 2 



f K/fc) = 0. 



v — U) r / k 



(60) 



(61) 



The second condition (jBlj) imposes that the phase velocity 
ui r /k = v ex t is equal to a velocity where f(v) is extrcmum 
(f'( v ext) = 0). The first condition (|60p then determines 
the wavenumbcr(s) k c corresponding to marginal stability. 
It can be written 



AttG [ + °° f'(v) , 
1 + ^^ / — dv = 0, 



k 2 . 



v - v ext 



(62) 



where the principal value is not needed anymore. The 
wavenumber(s) corresponding to marginal stability are 
therefore given by 

1/2 



-4ttG 



/» 
v - v ext 



dv 



(63) 



13 There exists less simple cases were explicit solutions can be 
obtained. We should mention in this respect the extensive set of 
closed-form solutions for generalized Lorentzian distributions 
due to Summers & Thorne (1991) and Thorne & Summers 
(1991). They correspond to Tsallis (polytropic) distributions 
of the form (|158|l with negative index n < — 1. 



Finally, the pulsation(s) corresponding to marginal stabil- 
ity are u) r = v ex tk c and we have Sf ~ e _8Wr *. Note that 
the distribution f(v) can be relatively arbitrary. There 
can be pure oscillations (cj = uj t ^ 0) only if f(v) has 
some extrema at v ^ 0. If f(v) has a single maximum at 
v = 0, then u r = (implying w = 0) and the condition of 
marginal stability becomes 



-4ttG 



dv 



1/2 



(64) 



3.5. The Nyquist method 



To determine whether the distribution f — f{v) is stable 
or unstable for a perturbation with wavenumber k, one 
possibility is to solve the dispersion relation (|52|) and de- 
termine the sign of the imaginary part of the complex pul- 
sation. This can be done analytically in some simple cases 
(see Sees. 13.81 0] and [5]). We can also apply the Nyquist 
method introduced in plasma physics. This is a graphical 
method that does not require to solve the dispersion rela- 
tion. The details of the method are explained by Nicholson 
(1992) and we just recall how it works in practice. In 
the e-plane, taking u>i = 0, we plot the Nyquist curv 
(e r (fc, oj r ), 6i(k, u) r )) parameterized by uj r going from —00 
to +00 (for a given wavenumber k). This curve is closed 
and always rotates in the counterclockwise sense. If the 
Nyquist curve does not encircle the origin, the system 
is stable (for the corresponding wavenumber k). If the 
Nyquist curve encircles the origin one or more times, the 
system is unstable. The number TV of tours around the ori- 
gin gives the number of zeros of e(fc,cj) in the upper half 
w-plane, i.e. the number of unstable modes with u>i > 0. 
The Nyquist method by itself docs not give the growth 
rate of the instability. 

Let us consider the asymptotic behavior of 
(e r (k, cj r ),€i(k, u r )) defined by Eqs. and 
for u> r ±00. Since f(v) is positive and tends to zero for 
v — > ±00, we conclude that ti(k, u) r ) — > for uj r — ¥ ±00 
and that 6i(k,u> r ) > for ui r — > —00 while 6i(k,ui r ) < 
for ui r — > +00. On the other hand, integrating by parts in 
Eq. (|5"8j) . we obtain 



4ttG r+oc 
e r (k,u} r ) = 1 + -p-" 



/(") 



(v — ui r /k) 



(65) 



provided that /(«) decreases sufficiently rapidly. 
Therefore, for uj r — > ±00, we obtain at leading or- 
der 



,, , 1 AirGp 
e r (fc,w r ) ~ H 5—, 



(uj r — > zboo) 



(66) 



In particular, e r (k,uj r ) — > 1 + for co r — > ±00. From these 
results, we conclude that the behavior of the Nyquist curve 
close to the limit point (1,0) is like the one represented 
in Fig. [3] In addition, according to Eq. (|59|) . the Nyquist 



This curve is also called an hodograph. 
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curve crosses the a;-axis at each value oiu> r /k correspond- 
ing to an extremum of f(v). For u r /k = v ext , where 
v ext is a velocity at which the distribution is extremum 
(f'( v ext) = 0), the imaginary part of the dielectric func- 
tion €i{k,kv ex t) = and the real part of the dielectric 
function 



€ r {k, kv ex t) 



1 



4ttG 



v - v ext 



dv. 



(67) 



Subtracting the value f'(v ex t) = in the numerator of the 
integrand, and integrating by parts, we obtain 



e r (k,kv ext ) = 1 



4ttG 



fjVext) ~ f(v) 

(v - v ext ) 2 



dv. 



If VMax denotes the velocity corresponding to the global 
maximum of the distribution, we clearly have 



e r (k,kv Max ) = 1 



4ttG 



fjvMax) ~ f(v) 
(V ~ V M ax) 2 



dv < 1. 



(69) 



Furthermore, e r (k, kvMax) < for sufficiently small k. 
Therefore, by tuning k appropriately, we can always make 
the Nyquist curve encircle the origin. We conclude that 
a spatially homogeneous self-gravitating system is always 
unstable to some wavelengths. 

3.6. Single- humped distributions 

Let us assume that the distribution f(v) has a single max- 
imum at v = vq (so that f'(vo) = 0) and tends to zero 
for v — > ±oo. Then, the Nyquist curve cuts the ir-axis 
(e,(fc,u; r ) vanishes) at the limit point (1,0) where uj r — > 
±oo and at the point (e r (k, fci>o), 0) where u r /k = vq. Due 
to its behavior close to the limit point (1,0), the fact that 
it rotates in the counterclockwise sense, and the property 
that e r (k i kvo) < 1 according to Eq. ([6"9")h the Nyquist 
curve must necessarily behave like in Fig. [3] Therefore, the 
Nyquist curve starts on the real axis at e r (k,uj r ) = 1 for 
uj r — > — oo, then going in counterclockwise sense it crosses 
the real axis at the point e r (k, kvo) < 1 and returns on the 
real axis at e r (fc,u; r ) = 1 for uj r — > +oo. According to the 
Nyquist criterion exposed in Sec. 13.51 we conclude that a 
single-humped distribution is linearly stable with respect 
to a perturbation with wavenumber k if 



e r (k, kv ) 



1 



4ttG 



m 

v~v 



■dv > 0, 



(70) 



and linearly unstable if e r (k, kvo) < 0. The equality corre- 
sponds to the marginal stability condition ([62]) . Therefore, 
the system is stable iff 



/ J-^-dv 



1/2 



(71) 

-oo V-V 

where k c is the critical Jeans wavenumber for a stellar 
system. Note that an infinite homogeneous stellar system 
whose DF has a single hump is always unstable to suffi- 
ciently small wavenumbers. For the unstable wavenumbers 



k < fe c , there is only one mode of instability Wi > since 
the Nyquist curve rotates only one time around the origin. 
This stability criterion is valid for any single-humped dis- 
tribution. In particular, a symmetric distribution / = f(y) 
with a single maximum at Vq = is linearly dynamically 
stable to a perturbation with wavenumber k iff 



k > k c 



-4ttG 



dv 



1/2 



(72) 



We shall make the connection between the stability of an 
infinite homogeneous stellar system and the stability of 
the corresponding barotropic gas in Sec. 13.91 In particular, 
using identity we will show that Eq. (fT2"j) is equivalent 
to Eq. JTT 



3. 7. Double-humped distributions 

Let us consider a double-humped distribution with a 
global maximum at VMaxi a minimum at v m i n and a lo- 
cal maximum at v max . We assume that VMax < Vmin < 
Vmax- The Nyquist curve will cut the x-axis at the limit 
point (1,0) and at three other points (e r (k, kvMax), 0), 
(e r (k, kv m in), 0) and (e r (fc, kv max ), 0). We also know that 
the Nyquist curve can only rotate in the counterclockwise 
sense and that e r (k,kvMax) < 1 according to Eq. ([69]). 
Then, we can convince ourselves, by making drawings, of 
the following results, i fi^l 

(+ + +): e r (vMax) > 0, e r ("mm) > 0, e r ( 

Umax )>o, 

(H ): £ r {vMax) > 0, e r (v mln ) < 0, e r (v max ) < 0, 

( £r{vMax) < 0, £r( u mm ) < 0, e r (v max ) > 0, 

(H h) : £r(vMax) > 0, e r (v m i n ) < 0, e r ( )>0, 

the Nyquist curve does not encircle the origin so the sys- 
tem is stable. If 

( ): e r (v M ax) < 0, e r (v mm ) < 0, e r (v max ) < 0, 

(- + +): e r (vMax) < 0, e r ( ) > 0, e r (v max ) > 0, 

(+ + -): £r(vMax) > 0, £ r ( ) > 0, e r (v 

max ) < o, 

the Nyquist curve rotates one time around the origin so 
that there is one mode of instability. Finally, if 

) < 0, e r (v min ) > 0, e r (v 

max ) < o, 

the Nyquist curve rotates two times around the origin so 
that there are two modes of instability. Cases (+ + +), 

( ), ( — h +) and ( — I — ) are observed in Sec. [7] for 

an asymmetric double-humped distribution made of two 
Maxwellians. The other cases cannot be obtained from 
this distribution but they may be obtained from other 
distributions. 

If the double-humped distribution is symmetric with 
respect to the origin with two maxima at ±v* and a min- 
imum at v — 0, only three cases can arise. If 

(+ + +): 6r(v*) > 0, e r (0) > 0, 

(+-+): er(«*) >0, e r (0)<0, 
the Nyquist curve does not encircle the origin so the sys- 
tem is stable. If 

( ):e r «)<0, e r (0)<0, 



15 In the following, in order to simplify the notations, we note 
e r {vMax) for £ r (fc, kVMax) etc. 
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the Nyquist curve rotates one time around the origin so 
that there is one mode of instability. Finally, if 

(- + -): e P («.) < 0, e r (0) > 0, 
the Nyquist curve rotates two times around the origin so 

that there are two modes of instability. Cases (+++), ( 

— ) and ( — I — ) are observed in Sec. 13.11 for a symmetric 
doublc-humpcd distribution made of two Maxwcllians. 



3.8. Particular solutions of e(k,u) = 

We can look for a solution of the dispersion relation 
e(fc, w) = in the form uj = iuji corresponding to uj r = 0. 
In that case, the perturbation grows (uji > 0) or decays 
(uji < 0) without oscillating. For un > 0, the equation 
e(fc, uj) = becomes 



(73) 



Multiplying the numerator by v + iuji/k and separating 
real and imaginary parts, we obtain 



k' 2 



oo V + 



/'(«) 



dv = 0. 



oo V + 



(74) 



(75) 



If wc consider distribution functions f(v) that are sym- 
metric with respect to v = 0, Eq. (|75|) is always satisfied. 
Then, the growth rate u>i > is given by Eq. ([74"]). 
For uji < 0, the equation e(k,uj) = becomes 

AnG f +co f(v) .8tt 2 G ( iuA 

1 + — L dv + l — f [tj = °- (76) 

Multiplying the numerator by v+iuji/k and assuming that 
f(v) is even, we obtain 

1 + ™ r™ * + -0, (77) 



k 2 



OO V + 



which determines the damping rate uji < 0. 

Let us introduce the function K(x) denned by 

1 / ,+0 ° u/'ful 
» «/'(«) 



where 



/ I ./-oo ^ + xZ 



dv + 2mf(ix) \ {x < 0) (79) 



This function is normalized such that K(0) = 1. The dis- 
persion relations (|74l) and (1771) can then be written 



(81) 



where fc c is the maginal wavenumber corresponding to cj = 
given by Eq. (fS4"|) . The pulsation Wj(fc) is given by 



k, = k K u 



(82) 



Setting u = uJi/k, it can also be written in parametric 
form as 



(83) 



where u goes from — oo to +oo. 

Let us obtain some asymptotic expansions of these re- 
lations (valid for symmetric distributions): 

(i) Let us first consider the case Wj > and k — > 
corresponding to instability. Integrating Eq. ([74")) by parts, 
we obtain 



(84) 



Expanding the integrand in powers of kv/uji <C 1, we find 
that 



w 2 = 4ttGp - 3Tfc 2 - 



(fc->0), 



(85) 



with T = (u 2 ) (where we recall that v = v z in the present 
case). This expression can be compared with the corre- 
sponding expression (|B.16|) valid for a gas. This identifica- 
tion yields c 2 = 3T so that large wavelength perturbations 
in a collisionless stellar system correspond to one dimen- 
sional isentropic perturbations with index 7 = 3 in a gas 
(see Appendix iBj) . 

(ii) The case uji < and k — > +00 corresponding to 
stability cannot be treated at a general level because the 
result depends on the behavior of the distribution function 
for large velocities. The case of a Maxwellian distribution 
is specifically considered in Sec. S) 

(iii) Let us finally consider the behavior of the disper- 
sion relation ((74)) or (|77j) close to the point of marginal sta- 
bility k = k c . For u>i = 0, we obtain the critical wavenum- 
ber k c given by Eq. (|M|. Let us consider w, —> + and 
k —> k~ in Eq. ([74]) . We introduce the function 



F(x) 



+ OO 



vf(v) 



dv, 



(86) 



for any real x. For x — > 0, we have the Taylor expansion 
F(x) = F(0) + F'(0)x + ... with 



F(0) 



+00 11 



f'(v) 



dv, 



(87) 



(80) F'(x) = - 



2xvf'(v) 
(v 2 + x 2 ) 2 



dv = — 



+00 



xf"(v) 



dv, 



where we have used an integration by parts to get the 
last expression. Under this form, we cannot take the limit 
x — ^ in the integral because the integral is not convergent 
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for x = 0. However, if we write Eq. 
form 

r+oo fl{v) _ r(Q) 



F'(x) 



we obtain 



F'(0) = - lim xf"(0) 

x->0 



dv 



in the equivalent 

+ °° xf"(0) A 
Z2Z — 2 dv > 

(89) 



dv 



= -7r/"(0)sign(x). 



(90) 

Regrouping the previous results, we find that the disper- 
sion relation (|74"]l becomes for — ^ + : 



1 



AttG 



/'(«) 



T/"(0)^ 



0. 



(91) 



Similarly, the dispersion relation (|77|) becomes for Wj — > 
0": 



1 + 



4ttG 



/'(«) 



*/"(0)^ 



which takes the same form as Eq. (|9"TT) . Therefore, Eq. (f9"TT) 
is valid for — > whatever its sign. To leading order, we 
obtain^ 



1 - 



k 2 



(k -> fc c ). 



(93) 



1 4tt 2 G/"(0) V fc?, 

This formula leads to the following result. First of all, we 
recall that the mode of marginal stability that we consider 
corresponds to uj r /k = v ex t = 0, i.e. to the extremum 
value of the distribution function f(y) at v ext = 0. If the 
distribution is maximum at v = 0, so that /"(0) < 0, 
we find that the mode uj = iuJi is stable for k > k c and 
unstable for k < k c . Alternatively, if the distribution is 
minimum at v = 0, so that f"(Q) > 0, we find that the 
mode uj = iuii is stable for k < k c and unstable for k > k c . 
This result will be illustrated in connection to Fig. [TJJ] for 
the symmetric double humped distribution. 

In this section, we have discussed particular solutions 
of the dispersion relation of the form uj = iuji. Of course, 
there may exist other solutions to the equation e(k, uj) = 
with u r 7^ C FH . However, for single-humped distributions 
and unstable wavenumbers k < k c , the Nyquist curve en- 
circles the origin only once (see Sec. 13. 6[) so that, when it 
exists, the solution uj = iu>i with cj, > is the only so- 
lution of e(k,co) = (for singlc-humpcd distributions and 
stable wavenumbers k > k c , there may exist other solu- 
tions of e(k, uj) = with oj r ^ and Wj < as discussed 
in Sec. I4.4|) . Explicit solutions of the dispersion relation 
e(fc, uj) = with uj = iuji arc given in Sees. [4] and [5] for 
isothermal and polytropic distributions. 



16 In the derivation, we have assumed that /"(0) 7^ 0. If 
/"(0) = 0, we need to develop the Taylor expansion to the 
next order. 

17 The modes w = ui r + iuji with Wj > and ui r 7^ are 
sometimes called overstable. 



3.9. Equivalence between the Jeans instability criterion 
for a stellar system and the Jeans instability 
criterion for the corresponding barotropic gas 

Lynden-Bell (1962) first observed that the critical Jeans 
length for a stellar system described by a Maxwellian dis- 
tribution function is equal to the critical Jeans length for 
an isothermal gas if we replace the velocity of sound by 
the velocity dispersion in one direction. In this section, we 
provide the appropriate generalization of this result for an 
arbitrary distribution of the form / = f(v 2 ) with /' < 0. 

As indicated in Sec. 13.21 to any stellar system with 
a distribution function / = /(e) and /'(e) < 0, we can 
associate a corresponding barotropic gas with an equation 
of state p = p(p). According to Eq. ©, the condition of 
hydrostatic equilibrium can be written 



P'(p) 



(94) 



(92) cf(r) 



Now, using Eq. ((?]) and the fact that p = J /(e) dv, we get 
p(r) p(r) p(r) 



Jf'{e)dv j Ml 



-dv 



J — OO V z A 



(95) 



where g(r, v z ) = J /(e) dv x dv y . For a spatially homoge- 
neous system, we obtain the identity 



(96) 



+00 f'(y) 



dv 



where we have noted v for v z and / for g like in Sec. 13.31 
This identity is explicitly checked in Appendix [C] for the 
isothermal and polytropic distributions. Since f(v) is sym- 
metric with respect to v — and has a single maximum at 
v = 0, the Jeans instability criterion can be written (see 
Sec. 



r < k 



-AttG 



f'(v) 



dv. 



Using identity (j^B]) . it can be rewritten 



k 2 <ki = 



AnGp 



(97) 



(98) 



Therefore, the criterion of dynamical stability for a spa- 
tially homogeneous stellar system coincides with the 
criterion of dynamical stability for the corresponding 
barotropic gas (see Sec. 12. 2[) . We insist on the fact that 
this equivalence is valid for an arbitrary distribution func- 
tion of the form / = f(v 2 ) with /' < 0, not only for the 
Maxwellian. We conclude that: an infinite homogeneous 
stellar system with / = f(v 2 ) and /' < is dynamically 
stable with respect to a perturbation with wavenumber 
k if and only if the corresponding barotropic gas is dy- 
namically stable with respect to this perturbation. This is 
the proper formulation of the Antonov first law for spa- 
tially homogeneous systems: in the present case, we have 
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equivalence^. Of course, although the thresholds of stabil- 
ity/instability coincide, the evolution of the perturbation 
is different in a stellar system and in a fluid system (see 
Sees. H] and [5]). We should also emphasize that, in general, 
the velocity of sound c s in the corresponding barotropic 
gas is not equal to the velocity dispersion (w 2 ) 1 / 2 in the 
stellar system, except when f(v) is the Maxwcllian distri- 
bution leading to Lyndon-Bell's result. 

4. Isothermal stellar systems 

4.1. The equation of state 

We consider an isothermal stellar system described by the 
distribution function 



(99) 



where (3 = 1 /T is a pseudo inverse temperature. We justify 
here this distribution as a particular steady state of the 
Vlasov equation^!. The associated pseudo entropy is 



S[f] 



/ln(///„)drdv, 



(100) 



where fo is a constant introduced to make the term 
in the logarithm dimcnsionlcss (it will play no role in 
the following since it appears in an additional constant 
term). The distribution (|99|) is obtained by extremizing 
the pseudo entropy (|100[) at fixed mass and energy, writing 
5S — 05E — aSM = 0. For an isothermal distribution, the 
kinetic temperature (velocity dispersion) defined by Eq. 
(|48p is spatially uniform, and it coincides with the pseudo 
temperature: T(r) = T. The barotropic gas corresponding 
to the isothermal stellar system defined by Eq. (|M)l is the 
isothermal gas with an equation of state p(r) = p(r)T. The 
velocity of sound is also spatially uniform and it coincides 
with the velocity dispersion in one direction: c 2 (r) = T. 
The density is related to the gravitational potential by Eq. 
(|2"0)) . It can be obtained by integrating Eq. (f9T)l) on the ve- 
locities leading to Eq. ([20]) with A' = A(2w/ /3) 3 / 2 . It can 
also be obtained by using Eq. © with p(r) = p(r)T, or 
by extremizing the functional (|21[) at fixed mass (Chavanis 
2006a). Finally, combining Eqs (j9"9")l and (|2"0")) , we can ex- 
press the distribution function in terms of the density pro- 
file according to 



/ 



3/2 



p(r) e 



(101) 



18 This "equivalence" for the dynamical stability of a homo- 
geneous stellar system and the corresponding barotropic gas 
differs from the case of inhomogeneous systems where the lim- 
its of stability of stars and galaxies are generically different 
(see, e.g., Chavanis 2006a). 

19 The statistical equilibrium state of a self-gravitating sys- 
tem (resulting from a "collisional" relaxation) is also described 
by an isothermal distribution of the form (|99|) . In that case, 
/ = Ae~ Pme is the Boltzmann distribution, /3 = l/(k B T) 
is the inverse thermodynamic temperature and Sb[/] = 
— ^ B J m m m^ r< ^ v * s ^ ne Boltzmann entropy (Padmanabhan 
1990). m 



Remark: the stability of box-confined isothermal stel- 
lar systems has been studied by Antonov (1962), Lynden- 
Bell & Wood (1968), Padmanabhan (1990) and Chavanis 
(2006a). 

4.2. The dispersion relation 

Let us now consider an infinite homogeneous isothermal 
system described by the Maxwcllian distribution function 
poll) with uniform density p(r) = p. The reduced distri- 
bution (ISTl) is 



/ 



1/2 



P e 



(102) 



The Maxwellian distribution has a single maximum at 
v = 0. Therefore, the condition of marginal stability 
dGTJi implies ui r = 0. From Eq. (j60|) . we find that the 
Maxwellian distribution is marginally stable for k = k c 
where we have introduced the critical wavenumber 



2 47rGp 



ki = 



T 



(103) 



According to the criterion (fTTj). the Maxwell distribution 
is linearly dynamically stable if k > k c and linearly dy- 
namically unstable if k < k c . The critical Jeans wavenum- 
ber (|103|) for an isothermal stellar system is the same as 
the critical Jeans wavenumber (|23[) for an isothermal gas. 
This is to be expected on account of the general result of 
Sec.EU 

The dielectric function (15^1) associated to the 
Maxwcllian distribution is 



e(k,u J ) = l- — l-) p 



3v 



■e-l 3 ^ dv. (104) 



Introducing the critical Jeans wavenumber (|103[) . it can 
be rewritten 



e (M) = i- 5? W( — 



where 



W{z) = 



/2tt 



2 dx, 



(105) 



(106) 



is the M^-function of plasma physics (Ichimaru 1973). We 
note that W(0) = 1. For any complex number z, we have 
the analytical expression 



W{z) = 1 - ze~^ 



f 

Jo 



e 2 dx 



(107) 



4.3. Growth rate and damping rate 

We look for particular solutions of the dispersion relation 
e(k,uj) = in the form uj = iuji where u>i is real. Using 
Eq. (fT0"7) , wc note that 



e(k,ioJi) 



k 2 



(108) 
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where we have introduced the function 



H{x) = W(ix) = 1 



— xe 2 erfc 



V2j 



(109) 



This function has the following asymptotic behaviors: (i) 
For x 0, H(x) = 1 - s/\x + .... (ii) For x +oo, 
H(x) = £(1 - + ...). (hi) For x -+ -oo, - 

—V27rxe~ (see Fig. [T]). Using e(k,iuji) = 0, the relation 
between and fc (for fixed T) can be written 



k 2 



UJ; 



fc 



0. 



or more explicit!; 



20 



fc 2 



— — — e^crfc \ — — 
2 k V V 2 k 



0. 



(110) 



(111) 



This equation can also be obtained directly from Eqs. (l53t 
and (|55p (see Appendix H)|) . If we set x = y/ftuji/k, we can 



rewrite Eq. (|110[) in the parametric form 
/77TT fc2 



fr(!c). 



(112) 



By varying x between — oo and +oo, we obtain the full 
curve giving w,; as a function of the wavenumber k (see 
Fig. [5]) • Since the time dependence of the perturbation is 



Sf 



the case of neutral stability uii = corresponds 



to k = fc c , the case of instability uii > corresponds to 
k < k c and the case of stability u>i < corresponds to 
k > k c . We have the asymptotic behaviors 



V4¥G^ ~ 2 fc 2 ' 




(fc/fc c -> 0), 



(fc/fc c 



(113) 



(114) 



(k/k c -> +oo). (115) 



Equation (|112l) provides a particular solution of the 
dispersion relation e(fc,w) = of the form uj = itu.i with 
w r = 0. The dispersion relation may have other solutions 
with uj r 7^ 0. However, for single humped distributions, we 
know that there is only one unstable mode with uji > 
for given k < k c (see Sec. l3.6p . Since the solutions uj = iuji 
given by Eq. (|112[) exist for any k < k c , we conclude that 
they are the only solutions in that range of wavenum- 
bers. Therefore, for the unstable wavenumbers k < k c , 
the perturbation grows exponentially rapidly without os- 
cillating. In other words, there are no overstable modes 
for the Maxwell distributional. This is the same behavior 

20 This relation is established by Binney & Tremaine (1987) 
for Ui > 0. The present analysis shows that it is also valid for 



tJi < 0. 
21 Binney 
method. 



Tremaine (1987) show this result by a different 



4- 



3 - 



x 



2 - 



2 4 
X 



10 



Fig. 1. The function H(x). 
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Fig. 2. Pulsation uj as a function of the wavenumber k 
for an isothermal stellar system. For k < kj, the system 
is unstable and uj = iuji with uii > 0. For k > kj, the 
system is stable. There exists many branches of solutions 
uj = uj r + iuii with uji < (see Sec. I4.4|) but we have only 
represented the branch corresponding to uj r = 0. We have 
also compared these results with the case of an isothermal 
gas. For k < fcj, the system is unstable and uj = ±iuji. For 
k > fc,/, the system is stable and uj — ±uj r . 



as in a fluid system (see Sec. [5J except that the growth 
rate uii > in the stellar system [see E q. 1 1 2 [) ] and in the 
fluid system [see Eq. (|l~7]) ] are differcntl 22 !. For the stable 
wavenumbers fc > fc c , the perturbations in a stellar sys- 
tem are damped exponentially rapidly (uji < 0). We have 
exhibited particular solutions (|112[) that arc damped with- 
out oscillating uj r = but these are not the only solutions 
of the dispersion relation. There also exists modes that are 
damped exponentially (uji < 0) and oscillate ui r ^ (see 
asymptotic results in Sec. 14. 4p . This form of damping for 
collisionless stellar systems is similar to the Landau damp- 



22 They only coincide for a cold system T = c 2 = (see Sec. 
or for fc = (see Sec. l3^)l . 
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ing for a plasmam. The situation is very different in a fluid 
system. In that case, the stable modes with wavenumbers 
k > k c correspond to gravity-modified sound waves that 
propagate without attenuation (cj r ^ 0, cjj = 0). 

The pulsation of the perturbations in an infinite ho- 
mogeneous isothermal stellar system is plotted in Fig. [5] 
as a function of the wavenumber k. For comparison, we 
have also indicated the pulsation of the perturbations in 
an infinite homogeneous isothermal gas. 

4.4. Other branches for k — > +oo 

Let us solve the dispersion relation e(k,u) = for an 
isothermal distribution in the limit k — > We look 

for a solution of the equation e(fc,w) = of the form 
u> = uj r + iixii with uii < (damping) and U)i ^> to r . This 
corresponds to heavily damped perturbations. We shall 
check this approximation a posteriori. Using Eq. (|55[) for 
uji < , Eq. can be written 



AttG 



+ oo 



/'(«) 



dv 



1tt 2 G 
k 2 



= 0. 



(116) 



If f(v) decreases like e @ v2 1 2 for real v — > ±oo, then for 
complex u> = u> r + i^i with uji 3> io r , f'{uj /k) will increase 
like e^/ 2fc . Therefore, to leading order, the foregoing 
equation reduces to 

8tt 2 G 



(117) 



Separating real and imaginary parts, we obtain two tran- 
scendant equations 



Re 



Im 



k 2 



-if 



= -I; 



'J' 



= 0, 



(118) 



(119) 



which crucially depend on the form of the distribution. 
For the Maxwcllian (|102j) . they can be rewritten to leading 
order in the limit Ldi/uj r 3> 1 as 

8tt 2 G ( 1/2 
2tt 



fc 3 



pfee& cos ( I = -1, 



k 2 



sin 



f3uj r uji 
k 2 



= 0. 



(121) 



Equation (|121l) implies (3u} r uji/k 2 = mix. Eq. (|120|) will 
have a solution provided that m is even. Then, Eq. (|120l) 
gives 

8tt 2 G / \ 1/2 ^? 



1. 



(122) 



In plasma physics, for the Maxwellian distribution, there 
is no solution to the dispersion relation of the form ui — iuii. 
The pulsation uj r is non zero (see Sec. 18.311 . 
24 We here adapt the method of plasma physics developed by 
Landau (1946), Jackson (1960) and Balescu (1963). 



which determines uii. By a graphical construction, it is 
easy to see that \u>i\ is an increasing function of k. For 
k — > +oo, we find the asymptotic behaviors 



VP 



7T 

-TO — 



1 



2 VP- 



(123) 



Since u>i/u) r ~ lnfc — > +oo, our basic assumption is sat- 
isfied. Therefore, for k > k c we have several branches of 
solutions parameterized by the even integer to. For to = 0, 
we recover the results of Sec. 14.31 

Remark: by analogy with plasma physics, we could also 
look for solutions of the dispersion relation ([5^]) of the 
form lo = ui r + iui with Ui <C ui r - This corresponds to 
weakly damped perturbations. In plasma physics, these 
solutions are valid for k — s- and lead to the usual Landau 
damping formula. However, a self-gravitating system is 
unstable for k < k c . Furthermore, it is easy to show that 
there is no solution of that form to Eq. (f512j) whatever the 
form of the distribution f(v) and the wavenumber k. This 
implies that for attractive interactions (like gravity) the 
perturbations are unstable for k < k c and heavily damped 
for k > k c while for repulsive interaction (like plasmas) 
they are weakly damped for k < krj and heavily damped 
for k > kjj. 

4.5. Nyquist curve 

It will be convenient in the following to work with di- 
mensionlcss quantities. We introduce the dimcnsionless 
wavenumber and the dimensionless pulsation 



AnGp 
Tk 2 ' 



O 



V^G~P 



(124) 



Noting that -^//Jw/fc = y^fJ, the dielectric function (|105[) 
can be rewritten 



e^o) = i- v w(Vnto)- 



(125) 



When Qi = 0, the real and imaginary parts of the di- 
electric function e(rj,Q r ) = e r (r/, tt r ) + iei{rj,VL r ) can be 
written 



e r (?7, fi r ) = 1 - ryW r (Vv^r) , 
ei(r), fl r ) = -rjWi (Vn^r) , 



with 

W r (z) 

Wi(z) 



dx, 



(126) 
(127) 

(128) 
(129) 



where z is here a real number. The condition of marginal 
stability corresponds to e r (?7, Sl r ) = ei(r],Q, r ) = 0. The 
condition ei(r?, fl r ) = 0, which is equivalent to f'(y/rjfl r ) = 
0, implies 57 r = 0. Then, the condition e r (r], f2 r ) = leads 
to t] = rj c with 



1. 



(130) 
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Fig. 3. Nyquist curve for the Maxwellian distribution 
(|102p . The DF is stable for perturbations with r\ < rj c 
(k > k c ), marginally stable for perturbations with 77 = n c 
(k = k c ) and unstable for perturbations with 77 > rj c 
(k < k c ). We have taken 77 = 2, 1,0.5 from the outer to 
the inner curve. 



To apply the Nyquist method, we need to plot the 
curve (e r (r7, r2 r ), 6^(77, J7 r )) in the e-plane. For f2 r — > ±00, 
this curve tends to the point (1,0) in the manner described 
in Sec. 13.61 On the other hand, for f2 r = 0, it crosses the 
x-axis at (e r (r],0) = 1 — 77, 0). The Nyquist curve is rep- 
resented in Fig. [3] for several values of the wavenumber 77. 
For 77 < 1 (i.e. k > k c ), the Nyquist curve does not encir- 
cle the origin so that the Maxwellian distribution is stable. 
For 77 > 1 (i.e. k < k c ) the Nyquist curve encircles the ori- 
gin so that the Maxwellian distribution is unstable. For 
77 = 1 (i.e. k = k c ) the Nyquist curve passes through the 
origin so that the Maxwellian distribution is marginally 
stable. The Nyquist method provides a nice graphical il- 
lustration of the Jeans instability criterion for an infinite 
homogeneous stellar system. 



5. Stellar polytropes 

5.1. The equation of state 

We consider a stellar polytrope (or polytropic galaxy) de- 
scribed by the distribution function 



/ = 



,% ~ 1) ' 
A* e 

q 



q-i 
+ 



(131) 



where /3 = 1/T is a pseudo inverse temperature and q is 
a parameter related to the traditional polytropic index n 
by 



1 



9-1 



(132) 



This relation is plotted in Fig. [Hand specific values consid- 
ered in the sequel are highlighted. We justify here the poly- 
tropic distribution function (|13ip as a particular steady 



state of the Vlasov equation^. The associated pseudo- 
entropy is 



S 



1 



9-1 



(133) 



where fo is a constant introduced for reasons of homogene- 
ity (it will play no role in the following since it appears in a 
term proportional to the mass that is conserved). The DF 
(|13ip is obtained by extrcmizing the pseudo entropy (|133|) 
at fixed mass and energy, writing 5S — f35E — aSAI = 0. 
The condition that C(f) must be convex imposes q > 0. 
On the other hand, we shall assume that /(e) is a decreas- 
ing function of the energy so that /? > 0. It is important to 
note that 1 //3 does not represent the kinetic temperature 
(or velocity dispersion) of the polytropic distribution (see 
Sec. EH} . 
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Fig. 4. The relation between the polytropic index 71 and 
the parameter q. Some particular values (q, n) are indi- 
cated for reference. 



25 Some authors (Plastino & Plastino (1993), Lima et al. 
(2002), Silva & Alcaniz (2004), Lima & de Souza (2005), 
Taruya & Sakagami (2003a), Leubner (2005), Kronberger et 
al. (2006), Du Julin (2006)) have interpreted the polytropic 
distribution (|131|) and the functional (|133|l in terms of Tsallis 
(1988) generalized thermodynamics. Here, we use a more con- 
ventional approach (Ipser 1974, Ipser & Horwitz 1979, Binney 
& Tremaine 1987, Chavanis 2006a). We interpret the poly- 
tropic distribution ()13ip as a particular steady state of the 
Vlasov equation and the functional (|133l) as a pseudo entropy. 
Its maximization at fixed mass and energy provides a con- 
dition of nonlinear dynamical stability with respect to the 
Vlasov equation (see Sec. 13. 1|) . not a condition of "general- 
ized thermodynamical stability". In particular, as argued by 
Chavanis & Sire (2005) and Campa et al. (2008), the insta- 
bilities reported by Taruya & Sakagami (2003b) correspond to 
Vlasov dynamical instabilities, not "generalized thermodynam- 
ical instabilities". In the present context, the analogies with 
Tsallis thermodynamics are purely coincidental. They are the 
mark of a thermodynamical analogy (Chavanis 2006a). Tsallis 
generalized thermodynamics applies in different contexts (see 
Chavanis 2008a). 
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Wc need to distinguish two cases depending on the sign 
of q — 1. For q > 1 (n > 3/2), the distribution function 
can be written 



f = A(e m -ey- 1 



+ > 



(134) 



where we have set A = [(3(q— l)/g]~ and e m = qp/[(3(q — 
1)]. Such distributions have a compact support since they 
vanish at e = e m . For e > e m , we set / = 0. Therefore, 
the notation [x]+ = x for x > and [x]+ = for x < 0. 
At a given position, the distribution function vanishes for 
v > v m (v) = y / 2(e m - $(r)). For q -> 1 (n -> +oo), we re- 
cover the isothermal distribution ([TO]) and for n = 3/2 the 
distribution function is a step function (see Sec. 15.51) . This 
is the distribution function of a Fermi gas at zero tempera- 
ture describing classical white dwarf stars (Chandrasckhar 
1942). For < q < 1, the distribution function can be 
written 



f = A(e Q + e) — , 



(135) 



where we have set A = [fi(l—q)/q]~ and e = qp/[/3(l — 
q)]. Such distributions are defined for all velocities. At a 
given position, the distribution function behaves, for large 
velocities, as / ~ u 2 /w -1 ) ^ w 2 ™ -3 . We shall only consider 
distribution functions for which the density p = f f dv 
and the pressure p =_§ J /i> 2 <iv are finite. This implies 
3/5 < g < 1 (n < -10. 

Let us now determine the equation of state of the 
barotropic gas corresponding to a stellar polytrope. For 
n > 3/2, the density and the pressure can be expressed as 

,r(3/2)r(n-l/2) 



4 7 rv / 2A(e m - $) r 



P 



-i-47r\/L4(e ro - $)' 



r(» + i) 

. x r(3/2)r(n-l/2) 



(136) 



(137) 



n + l - "' r(n + l) 

where T(x) denotes the Gamma function. For n < — 1, the 
density and the pressure can be expressed as 

r(-n)r(3/2) 



p = 4irV2A(e + 



T(3/2 - n) 



n + 1 1 (3/2 — n) 



(138) 



(139) 



Eliminating the gravitational potential between the ex- 
pressions (fl"56 > (fi"37 ]l and (fI55 ]) -([T59 ll . one finds that 



«--^H4$*f p»<-.). <-» 

Therefore, a stellar polytrope has the same equation of 
state (|140[) as a polytropic star. However, they do not have 
the same distribution function (compare Eq. (|131[) to Eq. 
(|4"|)) except for n — > oo corresponding to the isothermal 
case. The density is related to the gravitational potential 
by Eq. (f25|) . It can be obtained by integrating Eq. (|131[) on 
the velocity leading to Eqs. (|136[) and (|138[) from which, 
using Eqs. pH]) and (T4"2j) . we deduce Eq. {53} with A = 
e m /(K(n + 1)) for n > 3/2 and A = e /(-K(n + 1)) for 
n < —1. It can also be obtained by using Eq. (O with 
Eq. (|140[) or by extremizing the functional ([2l)f at fixed 
mass (Chavanis 2006a). Note that Eqs. ((25)) and (|26| are 
similar to Eqs. (|131[) and (| 133[) with 7 playing the role 
of the parameter q and K playing the role of the pseudo 
temperature T — 1//3. Polytropic distributions (including 
the isothermal one) are apparently the only distributions 
for which /(e) and p( ( f ) ) have a similar mathematical form. 

Remark: isolated stellar polytropes have a finite mass 
iff3/2<n<5 and they are dynamically stable (Binney & 
Trcmaine 1987). The stability of box-confined polytropes 
is studied by Taruya & Sakagami (2003a) in the context 
of generalized thermodynamics and by Chavanis (2006a) 
in the context of Vlasov dynamical stability. 

5.2. Other expressions of the distribution function 

We can write the distribution function of stellar poly- 
tropes (|131[) in different forms that all have their own 
interest. This will show that different notions of "tem- 
perature" exist for polytropic distribution^!^]. 

(i) Pseudo temperature T = 1//3: as indicated previ- 
ously, the form (|13ip of the polytropic distribution directly 
comes from the variational principle (|36[) when we write 
the pseudo entropy in the form (|133|) . Therefore, /? = 1/T 
is the Lagrange multiplier associated with the conserva- 
tion of energy. It is called "pseudo inverse temperature" . 
Note, however, that T = 1/(3 does not have the dimension 
of a temperature (squared velocity). 

(ii) Dimensional temperature 8 = 1/6: we can define a 
quantity that has the dimension of a temperature (squared 
velocity) by setting b = fi/qp. If we define furthermore 
/* = p l ^ q ^ 1 \ the polytropic distribution (|131[) can be 
rewritten 



P = Kp y , 7 =1 + -, 
n 



where 
K 



(140) 



n + 1 



r(n + 1) 



/ = /* 



U^2A T ^T {n : m \ in > 3/2) (141 



26 If we allow /3 to be negative, then it is possible to con- 
struct stellar polytropes with index 1/2 < n < 3/2 which are 
mathematically well-behaved (see Binney & Tremaine 1987). 
However, for those polytropes, the distribution function in- 
creases with the energy (and diverges at e = e m ) so they may 
not be physical. 



!-&(«-!)< 



(143) 



Using Eqs. (|136D and f)138[) . the relation between the den- 
sity and the gravitational potential can be written 



P = P* 



1 - b(q- 1)$ 



(144) 



27 We recall that, for collisionless stellar systems, the mass 
of the stars does not appear in the Vlasov equation and the 
different "temperatures" that we introduce have the dimension 
of a velocity squared. 
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with 

p* = 



2^[^V /2 r ^S r( "7, 1/2) (»>3/2), 



T(n+1) 



3 - 2n\ 3/2 r(3/2)r(- 



r(3/2-n) 



in < -1) 



(145) 



(146) 



(iii) Polytropic temperature K: eliminating the gravi- 
tational potential between Eqs. (|134[) and ()136|) . and be- 
tween Eqs. f)135[) and (|138[) . we can express the distribu- 
tion function in terms of the density according to 



/ 



p(r) 



l/n 



v 2 /2 



(n + l)K 



-3/2 



Z = 4^V2 r(3/ ^ r(n iN 1/2) [X(n + l)] 3 / 2 (n > 3/2)(148) 
T(n + 1) 



n < -1). (149) 



This is the polytropic counterpart of expression (|101[) for 
the isothermal distribution. The constant K plays a role 
similar to the temperature T in an isothermal distribution. 
In particular, it is uniform in a polytropic distribution as is 
the temperature in an isothermal system. For that reason, 
it is sometimes called a polytropic temperature. 

(iv) Kinetic temperature T(r): for a polytropic distri- 
bution, the kinetic temperature (velocity dispersion) de- 
fined by Eq. (|4"5]) is given by 



T(r) = Kpiv) 



7-1 



(150) 



For an inhomogeneous stellar polytrope, the kinetic tem- 
perature T(r) is position dependent and differs from the 
pseudo temperature T = 1//3. The velocity of sound is 
given by 



cj(r) = Kjpiry- 1 = 7 T(r) 



(151) 



It is also position dependent and differs from the velocity 
dispersion (they differ by a factor 7). Using Eq. (|150p . the 
distribution function (|147[) can be written 



/ = B n 



B n 



P(r) 



[27rT(r)] 3 / 2 



v 2 /2 



(n + l)T(r) 



n-3/2 



r(n + i) 



T(n- 1/2) (n+ I) 3 / 2 

r(3/2 - n) 
T(-n)[-(n+ l)]3/2 



(n > 3/2) 
(n < -1). 



(152) 



(153) 



(154) 



Note that for n > 3/2, the maximum velocity can be ex- 
pressed in terms of the kinetic temperature by 



v m (r) = y/2(n + l)T(r). 



(155) 



Using T(z + a)/T(z) ~ z a for z — > +00, we recover 
the isothermal distribution (|M| for n — > +00. On the 
other hand, from Eqs. (|150[) and (f2"5"j) . we immediately 
get T(r) = K(X - (7 - l)$(r)/if 7 ) so that 



7- 1 
VT= V$. 

7 



(156) 



This shows that, for a stellar polytrope, the kinetic tem- 
perature (velocity dispersion) is a linear function of the 
gravitational potential 2 ^!. The coefficient of proportion- 
ality is related to the polytropic index by (7 — l)/ 7 = 
l/(n + 1) = 2(q — l)/(5<7 — 3). This relation can also be 
obtained directly from Eq. (|134[) [or Eq. (|135[) ] noting that 



(147) / = A(e m - $) 



n-3/2 



1 - 



v 2 /2 



-1 n-3/2 



(157) 



and comparing with Eq. p52p . 
5.3. The dispersion relation 

Let us now consider an infinite homogeneous polytropic 
stellar system described by the polytropic distribution 
function (|152[) with uniform density p(r) = p and uniform 
kinetic temperature T(r) = T. From now on, T = 1//3 
will denote the kinetic temperature (|4"5)l , not the Lagrange 
multiplier appearing in Eq. (|13ip . The kinetic tempera- 
ture is uniform because we have assumed that the density 
is uniform. The reduced distribution function (1511) i 



12!) 



f(v) = B n 



V2^T 



2(n+ 1)T 



-1/2 



(158) 



where p is the density, T = 1//3 = (v 2 ) is the velocity 
dispersion in one direction and B n is a normalization con- 
stant given by 



B n — 
B n = 



r(n + l) 



r(? l + l/2)(n+l) 1 /2' 

r(l/2-n) 
r(-n)[-(n+l)]i/2' 



1 



n < — 1. 



(159) 
(160) 



The polytropic distribution has a single maximum at 
v = 0. Therefore, the condition of marginal stability (|61[) 
implies w r = 0. From Eq. (|60[) . we find that the polytropic 
distribution is marginally stable for k = k c where we have 
introduced the critical wavenumber 



AttGp 



(161) 



28 For any spherical stellar system with / = /(e), we have p = 

and p = £>($) so that the kinetic temperature (velocity 
dispersion) T = p/p is a function T — T(&) of the gravitational 
potential. For a polytropic distribution function, this relation 
turns out to be linear. 

29 If we justify the distribution function (|158|) from the 3D 
distribution function (I152|l integrated on v x and v y , then it is 
valid for n > 3/2 and n < — 1. However, as far as mathematics 
is concerned, the distribution (|158[) is normalizable and has a 
finite variance in the range of indices n > 1/2 and n < — 1. 
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According to the criterion (fTTj) . the polytropic distribution 
is linearly dynamically stable if k > k c and linearly dy- 
namically unstable if k < k c . The critical Jeans wavenum- 
ber (|16ip for a stellar polytrope is the same as the critical 
Jeans wavenumber (|3ip for a polytropic gas. This is to be 
expected on account of the general result of Sec. 13.91 It 
should be stressed that the quantity that appears in the 
critical wavenumber (|161l) is the velocity of sound c 2 = 7T 
in the corresponding barotropic gas, not the velocity dis- 
persion T. They coincide for the Maxwcllian distribution 
(7 = 1), but this is not general. 

The dielectric function (j5"2")l associated to the poly- 
tropic distribution is 



4ttG p ( 1 

e{k,uj) = 1 — \ n- - 



k 2 y/2^T 



1 



1 



2 J (n+l)T 

n-3/2 



2(n+l)T 



■ dv. 



v 



(162) 



Introducing the critical wavenumber (|161[) . we obtain 

',.2 



where 

W n {z) = 



1 B n ( 1 
/2tt n \ 2 



•[1- 



ira-3/2 



2(n+l)J 



(163) 

-dx, 
(164) 



is a generalization of the VF-function of plasma physics. 
We note that W n (0) = 1. For n — > +00, we recover the 
PF-function (|106p . Equation (|164p is therefore a general- 
ization of this function to the case of polytropic distribu- 
tions. 



5.4. Growth rate and damping rate 

We look for particular solutions of the dispersion relation 
e(fc,w) = in the form uj = iuji where u>i is real. First, we 
note that 



(165) 



where we have introduced the function H n (x) = W n (ix). 
For x > 0, we have 



H n {x) 



\_Bn 

/2n n 



+00 t 2 [l 



in-3/2 



2(n+l)J + 



t 2 



Using e(k,iu>i) = 0, the relation between u)i and k (for 
fixed T) can be written 



-dt, 
(166) 



and for x < 0, we have 

H n { X ) : 



1 B n ( 1 

7= I Tl — — 

2tt n \ 2 



+00 +2M _ t 2 ira-3/2 
1 I 1 2(n+l)J + 



-27TX 



1 



2(n+ 1) 



n-3/2 



k 2 



0. 



(168) 



For n — ^ +00, we recover Eq. (|110p . If we set x = \ffiu)i/k, 
we can rewrite Eq. (|168p in the parametric form 



'H n (x) 



7 



k 2 

-j^=H n {x). 



(169) 



By varying x between —00 and +00, we obtain the full 
curve giving u>i as a function of the wavenumber k. Since 
the time dependence of the perturbation is Sf ~ e Wit , the 
case of neutral stability w,; = corresponds to k = k c , the 
case of instability cj{ > corresponds to k < k c and the 
case of stability u)i < corresponds to k > k c . The dis- 
cussion of the different regimes is similar to the one given 
in Sec. 14.31 For k — > k c , using the general approximate 
formula (f9"3")l . the pulsation u>i is given by 
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(170) 
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Fig. 5. Pulsation ui as a function of the wavenumber k 
for a stellar polytrope with n > 1/2 (we have represented 
n = 0.5,1,1.5,2,3,5,10,100). For k < k c , the system is 
unstable and uj = iuji with oji > 0. For k > k c , the system 
is stable. There exists many branches of solutions uj = 
uj r + iuji with uji < but we have only represented the 
branch corresponding to uj r = 0. We have also compared 
these results to the case of a polytropic gas. For k < k c , 
the system is unstable and uj = ±iuji. For k > k c , the 
system is stable and uj = ±w r . The instability occurs for 



k < k 



(poly) 



k 



(iso) 



Apoly) . , (iso) 



i \fl where 7 = 1 4- 1/n. For n > 1/2, 



(167) 



The pulsation of the perturbation in an infinite homo- 
geneous stellar polytrope is plotted in Figs. [S] and [5] as a 
function of the wavenumber k. For comparison, we have 
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Fig. 6. Pulsation w as a function of the wavenumber fc 
for a stellar polytrope and a polytropic gas with n < — 1 
(we have represented n = —1.5, —2, —3, —5, —10, —100) 
The instability occurs for k < kc 
7 = 1 + 1/n. For n < 



-1 ki poLy ^ > ki iso * > 



= fcc JS °Vv / T' w here 



also indicated the pulsation of the perturbation in an infi- 
nite homogeneous polytropic gas. For n = oo, we recover 
the isothermal case shown in Fig. [5J 

5.5. Particular cases 

The case n = 3/2 deserves a particular attention. In that 
case, the distribution function /(r, v) is a step function so 
that / = 77o for v < v m (r) = v / 2(e m - $(r)) and / = 
otherwise. This corresponds to the distribution function 
of the self-gravitating Fermi gas at T = which de- 
scribes classical white dwarf stars (Chandrasekhar 1942). 
The density is p = ^TJoV^ and the pressure p = ||r/oW^. 
Eliminating v m from these two relations, we get a poly- 
tropic equation of state p = Kp 5 / 3 with n = 3/2, 7 = 5/3 
and K = ^(j^-) 2 ^ 3 - The kinetic temperature is T = ^v^ 
and the velocity of sound is c? a — For an infinite and 
homogeneous medium, the reduced distribution function 
P58p corresponding to n = 3/2 is a parabola 



(171) 



The critical Jeans wavenumber (f72j) is fc 2 = 12-KGp/bT = 
YlnGpIv^ which is fully consistent with the expression 
(fPoTj) with 7 = 1 + 1/n = 5/3. For n = 3/2, we can 
obtain an explicit expression of the function (|166[) - (|167[) . 
For x > 0, we have 



H 3/2 {x) = 1 



x f V5 
— ■= arctan — 
V5 \x i 



(172) 



and for x < 0, we have 



H 3/2 (x) = 1 - -j= 




The index n = 1/2 is also special and corresponds to 
the water-bag model. In that case, the reduced distribu- 
tion f(v) is a step function so that / = 770 for \v\ < v m and 
/ = otherwise. The amplitude 770 is determined by the 
density according to the relation p = 2rjQV m . The kinetic 
temperature is T — (v 2 ) = \v 2 n . The derivative of the 
distribution function is f'(v) = r/o[5(v + v rn ) — 6(v — v m )]. 
The critical Jeans wavenumber ((72^1 is fc 2 = AirGp/v 2 = 
AnGp/iT which is fully consistent with the expression 
(|16ip with 7 = 1 + 1/n = 3. For n — 1/2, we can ob- 
tain an explicit expression of the dielectric function (|162[) . 
We get 



with 

W l/2 = 



k 2 



kVT 



1 



3 Z 



(174) 



(175) 



The condition e(k,u>) = determines the pulsation. For 
k > fc c , the system is stable and the perturbation presents 
pure oscillations with pulsation 



±V3T(k 2 - fc 2 ) 1 / 2 . 



(176) 



For k < k c , the system is unstable and the perturbation 
has a growth rate (and a decay rate) given by 



uj = ±tV3T(k 2 c -k 2 )^ 2 . 



(177) 



We note that, for the water-bag distribution, the general 
asymptotic behavior (|55j) becomes exact for all k < k c . 
We also note that for the specific index n = 1/2 (7 = 3), 
the dispersion relation in a stellar system takes the same 
form as in a gas (see Sec. 0). 

5.6. The Nyquist curve 

Introducing the dimensionless wavenumber and dimen- 
sionlcss pulsation (|124[) . the dielectric function f| 163[) can 
be rewritten 



e(r,,fl) = l--r 1 W n ( y /?jQ). 

7 



(178) 



When f2i = 0, the real and imaginary parts of the di- 
electric function £(77, Q r ) = e r (r), fi r ) + ^(77, fi r ) are given 

by 



e r (r ] ,n r ) = l--riW< : n H^rjn r ). 
7 

7 



(179) 
(180) 



with 



V27T 77 V 2 



(173) xP 



+°° x[l 
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-dx, 



(181) 
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1 - 



2(n + l) 



1-3/2 



(182) 



where z is here a real number. The condition of marginal 
stability corresponds to e r (r/, fl r ) = e,; (77, fi r ) = 0. The 
condition £4(77, f2 r ) = 0, which is equivalent to f'(^/rjQ r ) = 
0, implies Q r = 0. Then, the relation e r (?7, 51 r ) = leads 
to ?/ = r\ c with 



% = 7- 



(183) 



To apply the Nyquist method, we need to plot the curve 
(e r (?7, f2 r ), £4(77, S^r)) in the e-plane. We have to distinguish 
different cases according to the value of the index 71. A 
general discussion has been given by Chavanis & Dclfini 
(2009) in the context of the HMF model. This discussion 
can be immediately transposed to the present context. 




Fig. 7. Symmetric doublc-humpcd distribution made of 
two Maxwellians with separation y (for a given tempera- 
ture T). If y > 1, the distribution has two maxima at ±V* 
and one minimum at Vq = while for y < 1, it has only 
one maximum at Vq = 0. 



6. The symmetric double-humped distribution 

6.1. Determination of the extrema 

We consider a reduced distribution function (jSTj) of the 
form 




(184) 



This is a symmetric doublc-humpcd distribution corre- 
sponding to the superposition of two Maxwellian distri- 
butions with temperature T = 1//3 centered in v a and 
— v a respectively (see Fig. [7]). This distribution models 
two streams of particles in opposite direction. The av- 
erage velocity is (v) = and the kinetic temperature 
Tf-in = (v 2 ) = T + v 2 . The velocitie(s) vq at which the 
distribution function f(v) is extremum satisfy f'(vo) = 0. 
They are determined by the equation 



-2f)v a v _ V °- ~ ^0 
V a + Vq ' 



(185) 



We note that vq G [— v a , +v a ]. Introducing the dimension- 
less velocity and the dimensionlcss separation 



V = y/fiv, V a = ^fpV a , 

Eq. (|185p can be rewritten 
1 , fV a + Vo 



1 



hi 



2V a V \Va-V . 
It is convenient to introduce the variables 

* a 



(186) 



(187) 



(188) 



For a fixed temperature T, x plays the role of the velocity 
vq at which the distribution is extremum and y plays the 
role of separation v a . Then, we have to study the function 



(189) 



for x e] — 1, +1[. This function is plotted in Fig. [S] It has 
the following properties: 



y(-x) = y(x), 



y(o) = 1. 



(190) 
(191) 
(192) 



The extrema of the distribution function (|184p can be de- 
duced from the study of this function. First, considering 
Eq. (|187|) . we note that f(v) always has an extremum at 
vq = 0, for any value of (3 and v a . This is a "degenerate" 
solution of Eq. (|189|) corresponding to the vertical line 
x — in Fig. [5J On the other hand, if y > 1, i.e f3v 2 > 1, 
there exists two other extrema vq = ±u* where tj* = v a x* 
with x* = y^d^v 2 ). 

In conclusion, for a given temperature T: 

• if y > 1 (v 2 > T), the distribution function f(v) has 
two maxima at vq = ±u* and one minimum at vq = 0. 

• if V < 1 ( w a — ^)) t ne distribution function f(v) 
has only one maximum at vo = (the limit case /3V 2 . = 1 
corresponds to /"(0) = 0). 

6.2. The condition of marginal stability 



Introducing the dimensionless variables (|124[) and (|188|) . 
the dielectric function associated to the symmetric double- 
humped distribution (|184[) is 



(193) 



where W(z) is defined in Eq. (| 106|) . For a fixed temper- 
ature, r] plays the role of the wavenumber k, Q plays the 
role of the pulsation w and y plays the role of the separa- 
tion v a . When 57^ = 0, the real and imaginary parts of the 
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Fig. 8. The function y(x) for the symmetric double- Fig. 9. rj = l/W r {y/y) as a function of y. 
humped distribution. 



dielectric function e(r?, fi r ) = 
written 



e r (?7, f2 r ) + iei(r], fl r ) can be 



1.708, W r (y/y) is negative so the branch 
For y — > 0, we have 



(200) 



ei{r),n r ) 



[Wi(y/rjn r - y/y)) + Wi(^rjn r + y/y))} , 



where W r {z) and Wi(z) are defined in Eqs. ([T28 > 1(129 )1 
where z is here a real number. The condition of marginal 
stability corresponds to e r (r], f2 r ) = 6,(77, fi r ) = 0- The 
condition (77, f2 r ) = is equivalent to 



For y > y„ 

Vc (y) exists only for y G [0, y n 

Vv))], ^(o) = 1. 

(I94THS result is to be expected since, for v a = 0, the distri- 
bution (|184p reduces to the Maxwellian. We thus recover 
the critical wavenumber (|130[) . 
In conclusion: 

• if y < y m axi there exists a critical wavenum- 
ber r]c°\y) determined by Eq. Q198P corresponding to a 
marginal mode (f2 r = 0, Qi = 0). 

• if y > y m ax, there is no marginal mode (f2 r = 0, fij = 

0). 



(195) 



The condition e r (r], Q, 
V 



1 - \ [Wr(^r 



Vv)) 



leads to 



(196) 

v / y))] = °- ( 197 ) 



Therefore, according to Eq. f|196[) . the phase velocity 
y/rj£L r is equal to a velocity Vq at which the distribution 
(I184|) is extremum. The second equation (|197[) determines 
the value(s) rj c (y) of the wavenumber at which the distri- 
bution is marginally stable. 



6.2.1. The case u r 







Let us first consider the value O r = that is solution of 
Eq. (|196p for any v a and (3. In that case, Eq. (|197p becomes 

»< 0) (») 



where we have used W r {— x) = W r (x). For given separa- 
tion y, this equation determines the wavenumber rji°\y) 
corresponding to a mode of marginal stability with fl r = 
0. The function defined by Eq. (fl~98j) is plotted in Fig. [9] 
It diverges at y = y m ax = z 2 c where z c = 1.307 is the zero 
of W r (z) (see Appendix A of Chavanis & Delfini 2009). 
Then, using W^{z c ) = -l/z c , we find from Eq. (fH?g]) that 

max 



(jj -> ymax)- 



(199) 



The curve (y) corresponding to the marginal mode 
with zero pulsation fi r = is plotted in Fig. [TUJ 

6.2.2. The case uj r / 

We now consider the cases where y/rjfl r = ±V* is solution 
of Eq. (|196[) for y > 1. To determine the wavenumbcr(s) 
at which the distribution (|184[) is marginally stable, we 
have to solve 



1 - j [W r {V* - y/y)) + W r (V* + y/y))] = 0, 
where V* is given by 



1 



1 



y/V + V* 
2V.y/y-\y/y-V. 



In 



(201) 



(202) 



Eliminating V* between these two expressions, we ob- 
tain the critical wavenumber(s) rfc^\y) as a function of 
y. However, it is easier to proceed differently. Setting 
x = V*/V a = Vt/y/y, we obtain the equations 



y 



2x 



1 



[W r (y/y(x-l)) + W r (y/y(x + l))] 



(203) 
(204) 



For given x, we can obtain y from Eq. (|203p [see also Fig. 
[8] and rj from Eq. (|204p . Varying x in the interval ] — 1, 1[ 
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yields the full curve 77^ (y) . By symmetry, we can restrict 
ourselves to the interval x £ [0, 1[. 
For x = 0, we have y = 1 and 



# } (1)=^ = 



f 



W r {\) 



= 3.633. 



(205) 



The branch rji^iy) starts at the point (1, 77*), correspond- 
ing to y^rjflr = ±14 = (i.e. x = 0). This point is at 
the intersection between the branch r]i°\y) along which 
Q r = and the line y = 1 separating the regions where 
the distribution has one or two maxima. 

For x — > 1, we have y ~ — i ln(l — ar) — > +oo and 



1 

2^ 



(y -> +oo). 



(206) 



This result can be understood simply. For 77 — > +oo, the 
two humps are far away from each other so that the critical 
wavenumber 77^ coincides with the critical wavenumber 
of a single Maxwcllian since they do not "see" each other. 
Noting that the mass of a single hump is M/2, the corre- 
sponding critical wavenumber is k c = 4nG(p/2)/T leading 
to 4 ±] = 2. 

In conclusion: 

• if y > 1, there exists a single critical wavenumber 
Vc^\y) determined by Eqs. (|203[) - (|204p corresponding to 
a marginal mode (f2 r = ztVt/y/rf^, f2, = 0). Note that the 
modes fl r = +V*/y/ffc and Q, r = —V^/^/rj^ are degener- 
ate. This degeneracy can be raised by a small asymmetry 
(symmetry breaking) in the distribution (see Sec. [7]). 

The curve (y) corresponding to the marginal mode 
with pulsation il r = ±.V*/y/rfc is plotted in Fig. ITUl 



6.3. The stability diagram 

The critical wavenumbers rj c (y) corresponding to marginal 
stability determined previously are represented as a func- 
tion of the separation y in Fig. [TDJ We have also plotted 
the line y = 1. On the left of this line, the distribution 
has a single maximum at Vq = and on the right of this 
line, the distribution has two maxima at Vo = ±V* and a 
minimum at Vq = 0. In order to investigate the stability 
of the solutions in the different regions, we have used the 
Nyquist method. 

For y < 1, there exists one wavenumber ric \y) at 
which the DF is marginally stable. For r\ = T)^\y), the 
DF has one maximum at Vo = 0. The marginal pertur- 
bation does not propagate (O r = 0). By considering the 
Nyquist curves in this region (see Figs. I1H12[) , wc find 
that the DF is stable for 77 < i]i°\y) and unstable for 
r] > ??c 0) (jy). 

For y > y max , there exists one wavenumber (y) at 
which the DF is marginally stable. For 77 = r]c ^ {y), the DF 
has two maxima at Vq = ±V* and one minimum at Vo = 0. 
The marginal perturbation evolves with a pulsation O r = 
±V*/y/rfc. By considering the Nyquist curves in this region 



(- + -) 



Unstable 
(two modes) 




(x=lj 



Fig. 10. Stability diagram of the symmetric double- 
humped distribution (|184[) . The solid line (7y c )* corre- 
sponds to the critical line: below this line the DF is stable 
and above this line the DF is unstable. On the left panel 
(delimited by the dotted line), there is one mode of in- 
stability and on the right panel there are two modes of 
instability. The dashed line corresponds to y = 1: on the 
left of this line the DF has one maximum and on the right 
of this line the DF has two maxima and one minimum. 
The symbols in parenthesis like (x = 0) give the values 
of x or y that parametrize the marginal curves. The no- 
tations like ( — I — ) give the positions of the e r (v ext Ys in 
the Nyquist curves as defined in Sec. 13.71 
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Fig. 11. Nyquist curve for y < 1 and 77 < (specifically 
y = 0.5 and 77 = 0.5). The DF has only one maximum at 
Vq = 0. It is stable (with respect to this perturbation) 
because the Nyquist curve does not encircle the origin. 
Case (+). 



(see Figs. IT3"l[i"4"|) . we find that the DF is stable for 77 < 
Vc^\y) an d unstable for 77 > rjc (y). 

For 1 < y < y ma x, there exists two wavenumbers 
77^ (77) and (y) at which the DF is marginally stable. 
The DF has two maxima at Vo = ±V* and one minimum 
at Vo = 0. For 77 = rjc (y), the marginal perturbation does 
not propagate (f2 r = 0). For 77 = rjc (y), the marginal 
perturbation has a pulsation fl r = ±V*/y/rfc. By consid- 
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Fig. 12. Nyquist curve for y < 1 and ?y > T]c (specifically 
y = 0.5 and 77 = 6). The DF has only one maximum at 
Vo = 0. It is unstable (with respect to this perturbation) 
because the Nyquist curve encircles the origin. Case (-). 
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0.4 - 
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Fig. 13. Nyquist curve for y > y max and r) < r)c (specif- 
ically y = 2 and 77 = 1). The DF has two maxima at 
Vq = ±V* and one minimum at Vo =0. The DF is stable 
(with respect to this perturbation) because the Nyquist 
curve does not encircle the origin. Case (+ + +). 

ering the Nyquist curves in this region (see Fig. I15[) . we 
find that the DF is stable for 77 < r)c(y) and unstable for 

V > Vc^iv)- 

A few comments arc in order: 

1. In Fig.[T0l we explicitly see that the system is always 
unstable with respect to perturbations with sufficiently 
small k (the case of cold systems T = is treated in 
Appendix [EJ. This corroborates the general result given at 
the end of Sec. 13.51 More precisely, we see that the system 
is stable for perturbations with k > (k c )* (corresponding 
to the solid line) and unstable for perturbations with k < 
(fc c )*. Furthermore, we note that this critical wavenumber 
(fc c )* corresponds to a marginal perturbation where the 
phase velocity u/k coincides with the maximum of the 
velocity distribution: for v 2 a < T, this is vq = and for 
v 2 a > T, this is vq = ±u». 




-2 -1 1 



Fig. 14. Nyquist curve for y > y max and 77 > T)c (specif- 
ically y = 2 and rj = 6). The DF has two maxima at 
Vq = ±V* and one minimum at Vo — 0. The DF is unstable 
(with respect to this perturbation) because the Nyquist 
curve encircles the origin. Since it rotates twice around 
the origin, this implies that there are N = 2 unstable 
modes (u; r ,w,) with uji > 0. Case (- + -). 




e 



Fig. 15. Nyquist curve for 1 < y < y rnax and 77 > i]c 
(specifically y = 1.2 and r\ = 6). The DF has two maxima 
at Vo = ±14 and one minimum at Vb = 0. The DF is 
unstable (with respect to this perturbation) because the 
Nyquist curve encircles the origin. Since it rotates only 
once around the origin, this implies that there is N = 1 
unstable mode (w r ,u)i) with uji > 0. Case ( ). 



2. For 1 < y < y max and tjc < rj < rfc , there are 
two unstable modes since the Nyquist curve encircles the 
origin twice. One of these two modes is lo — iuii and the 
other is uj = w r + iu)i with u> r 7^ (overstable) . When we 
increase r\ and cross the marginal line rj = , the mode 
uj = iuji becomes stable according to the general result 
(|93[) . Indeed, for y > 1, the DF is minimum at v — 0. This 
is why there is only one unstable mode for rj > r]c . the 
mode uj = uj r + iuji with uj r ^ (overstable) . On the other 
hand, for y < 1, the mode uj = iuji is stable for r\ < r]i \ 
When we increase 77 and cross the marginal line 77 = rjc , 
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it becomes unstable according to the general result ([Ml) . 
Indeed, for y > 1, the DF is maximum at v = 0. 

3. In Fig. [TU1 we see that the critical Jeans length 
(A c )* associated to a double-humped DF (v a ^ 0) is al- 
ways larger than the critical Jeans length associated to 
the single- humped Maxwellian (v a = 0). This means that 
the presence of streaming in a purely stellar system has 
a stabilizing role. This is different if there exists a gas 
component in the system since the critical Jeans length is 
reduced by the relative motion of the gas and stars (Sweet 
1963). 

4. For a double-humped stellar system with v a < \/T, 
the critical Jeans length increases with v a due to the sta- 
bilization effect of the relative velocity. Alternatively, in 
a contra-streaming self-gravitating gas with v a < c s , the 
critical Jeans length decreases with v a and tends to zero 
as v a — > c s (Talwar & Kaha 1966, Ikeuchi et al. 1974). 
On the other hand, for a double-humped stellar system 
with v a > VT or for a contrastreaming self-gravitating 
gas with v a > c s , the critical Jeans length decreases with 
v a (but remains larger than the classical Jeans wavelength 
corresponding to v a = 0) and overstable modes appear. 

7. The asymmetric double-humped distribution 

7.1. Determination of the extrema 



We now assume that the reduced distribution (|51[) is an 
asymmetric double-humped distribution of the form 




r (»-»») 



(207) 



where T = 1//3 is the temperature of the Maxwellians and 
A is the asymmetry parameter (we assume here that A > 
1). This distribution is plotted in Fig. 1161 The symmetric 
case is recovered for A = 1. The average velocity is (v) = 
— [(A — 1)/(A + l)]v a and the kinetic temperature T^m = 
((« - (v)) 2 ) = T+ [4A/(A + l) 2 ]v 2 a . The velocities v at 
which the distribution function f(v) is extremum satisfy 
f'{vo) = 0. They are determined by the equation 



-2Pv a v = 1 Z 

A v a + v ' 



(208) 



We note that vq g] — v a , +v a [. Introducing the dimension- 
less velocity and the dimcnsionless separation (|186|) . Eq. 
(1208|) can be rewritten 



1 



1 



hi 



Va + V 



ln(A) 



2V a V \V a -V J 2V a V 
It is convenient to introduce the variables 

X = 7T> y = v a . 

v a 

Then, wc have to study the function 



y(x) = i ln (f 



ln(A) 

2x 



(209) 



(210) 



(211) 
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Fig. 16. Asymmetric double-humped distribution made of 
two Mawellians with separation y and asymmetry A > 1. 
If y > y*, the DF has one global maximum at Vq = V- < 
0, a minimum at Vq = V p > and a local maximum at 
Vo = V+ > 0. If y < y*, the DF has only one maximum at 
Vo = V- < 0. 



for x g] — 1, +1[. This function is plotted in Fig. [17] It has 
the following properties 



y(x) iln(l - x), (x-tl ), 



y(x)~--ln(l + x), (x^-l + ), 



y(x) 



ln(A) 

2x : 



(x -4 0). 



(212) 
(213) 
(214) 



Considering the negative velocities x < 0, we note that 
y > iff x < xo with 



1 



x 



1 



(215) 



Considering the positive velocities x > 0, we note that the 
curve y{x) is minimum at x = cc* where x* is solution of 



2x# 



l-x 2 



-In 



1 + x t 
l-x> 



ln(A). 



(216) 



This function is represented in Fig. [TH We note that Eq. 
(I216P has a unique solution x„ for each value of A > 1. 
Therefore, the function y(x) has a single minimum at x = 
x* . The value of this minimum is 



1 



1-xi 



> 1. 



(217) 



The extrema of the distribution /(V) can be determined 
from the study of the function (|211|) . If y < y*, the dis- 
tribution f(V) has a single maximum at Vq = V— < 0. 
If y > y*, the distribution f(V) has two maxima at 
Vq — V- < and Vo = V+ > and one minimum 
at Vq = V p > 0. These different values are given by 
Vo = V a y- 1 (V 2 ). 

In conclusion, for a given asymmetry A > 1 and tem- 
perature T: 
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Fig. 17. The function y(x) for the asymmetric doublc- 
humpcd distribution. 




Fig. 18. Evolution of x*, j/*, x c and y c as a function of A. 
The curves intersect each other at A* = 3.3105. 



• if y > j/*, the distribution function f(V) has a global 
maximum at Vq = V— < 0, a local maximum at Vb = 
V+ > and one minimum at Vo = V p > 0. 

• if y < y*, the distribution function /(V) has only 
one maximum at Vb = V— < 0. 



7.2. 77?e condition of marginal stability 

The dielectric function associated to the asymmetric 
double-humped distribution is 



e(fi) = 1 - 



V 



1 + A 



- y/y)) + zwyrjQ. + ^))\ , 



where W r (z) and Wj(z) are defined in Eqs. 1(125 - ([125 )1 
where z is here a real number. The condition of marginal 
stability corresponds to e r (r), fl r ) = 6^(77, Q r ) = 0. The 
condition ei(rj, Q r ) = is equivalent to 



/'(VW = o. 

The condition e r (rj, Cl r ) = leads to 



(221) 



1 - 



1 + A 



[W r (07fir - y/y)) + AW r (07f2 r + y/y))] = 0. 



(222) 



According to Eq. (|22 1|> . the phase velocity y/rj£l r is 
equal to a velocity Vq at which the distribution (|207p 
is extremum. The second equation (|222p determines the 
value(s) ?7c(y) of the wavenumber at which the distribu- 
tion is marginally stable. Therefore, we have to solve 



n 



1 + A 



[W r (V - y/y)) + AW r (V + y/y))] = 0, 



where Vq is given by 



1 



1 



In 



fVy + v °\ , ln (A) 



2y/Wo \y/V-V 



ly/yVo 



(223) 



(224) 



Eliminating Vq between these two expressions yields the 
critical wavenumber(s) T) c (y) as a function of y. However, 
it is easier to proceed differently. Setting x = Vq/ y/y, we 
obtain the equations 



y 



2x 



ln(A) 
2x ' 

+ A 



[W r (y/y(x - 1)) + AW r (^y(x + 1))] 



(225) 



(226) 



For given x, we can obtain y from Eq. (|225[) [see also 
Fig. \T7\ and r\ from Eq. (|226p . Varying x in the interval 
] — 1, 1[ yields the full curve r) c (y). We have three types of 
solutions. For x s] — l,xo], we obtain a branch rjc \y) 
where the pulsation of the marginal mode is negative: 
il r = V-/y/rj^ < corresponding to the global maxi- 
mum of f(v). For x €]0, x*], we obtain a branch (y) 
where the pulsation of the marginal mode is positive: 
Q r = V p / y/rfc > corresponding to the minimum of f(v). 

(+)/ 



where W(z) is defined in Eq. (|106f) . When Oj = 0, the real 
and imaginary parts of the dielectric function e(r?, Q r ) = 
e r (r], Q r ) + iei(j], Q r ) can be written 



e r (rj, fl r ) = 1 



n 



^[W r {y/r)£l r -y/y)) 

+AW r (v^o r + Vy))] J 



(218) For x G [x*, 1[, we obtain a branch t]c (y) where the pul- 
sation of the marginal mode is positive: Q r = V+/ y/ff c > 
corresponding to the local maximum of f(v). This leads 
to the curves reported in Figs. IT9"1 and I2U1 for two values of 
the asymmetry factor A. 

For x — > —1, y — > +oo and 



(219) v (,-) {y) 



1 



1 



ei{r], O r ) 



>7 



^[Wi(V^r--S^)) 

+Aw i (0?fir + Vy))] J 



(220) 



A 4A 2 y 
For x — > xo, y — > and 

7?<->fo)->l (y->0). 



(y -> +00). 



(227) 



(228) 
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This returns the critical wavenumber (|130[) associated 
with the Maxwellian distribution (|102[) corresponding to 
y = 0. For x — > 1, y — > +oo and 

4+\y) ~ 1 + A+ A(1 J A) + ■■■ (y^+oo). (229) 

Let x c and y c be determined by the equations 

ln(A) 



1 . fl + x, 
2x r VI- x 



W r (y/y2(x c 
For a; — > x c 

vi s) (y) « 



2x c ' 

1)) + AW r (Vlte(a: 
— > y c and 

> +oo (y -> j/ c ) 



1)) = 0. 



(230) 
(231) 

(232) 



where s = p if A < A* and s = + if A > A* as will 
become clear below. Note that there is no physical solu- 
tion to Eqs. (|225|l - (|226|l when < x < x c (rj would be 
negative) so that the branch r]c (y) starts at (y c ,+oo) 
corresponding to x = x c . The evolution of x c and y c with 
A is studied in Fig. Upland is compared to the evolution of 
x* and y*. The curves intersect each other at A* = 3.3105. 
For A < A*, x c < x* so that the stability diagram dis- 
plays three marginal branches r]c (y) , r]c (y) and r]c (y) 
(see Fig. \T^i . The branches r)c P \y) and ijc + \y) connect 
each other at (y„ , rj* ) corresponding to x = x* . At that 
point they touch the line y = y* separating distributions 
with one or two maxima. For A = A*, x c = x* so that the 
branch rji P \y) is rejected to infinity and only the branches 
rjc \y) and r)c (y) remain. For A > A*, x c > x* so that 
the phase diagram displays only two marginal branches 
7}i~\y) and ?] { c +) (y) (sec Fig. HQ). 



7.3. The stability diagram 

The critical wavenumbers rj c (y) corresponding to marginal 
stability determined previously are represented as a func- 
tion of the separation y in Figs. IT9l and l20l for two values 
of the asymmetry factor A. We recall that rjc (y) cor- 
responds to the wavenumber associated with a marginal 
mode with pulsation tt r = V- / ypffc < (global maximum 

of /), rji p \y) corresponds to the wavenumber associated 
with a marginal mode with pulsation £l r = V p / ' ^frf c > 

(minimum of /) and r]i + \y) corresponds to the wavenum- 
ber associated with a marginal mode with pulsation fi r = 
V+/y/Vc ^ (local maximum of /). We have also plotted 
the line y = y*. On the left of this line, the DF has a sin- 
gle maximum at V- < and on the right of this line, the 
DF has two maxima at V_ < and V+ > and a mini- 
mum at V p > 0. In order to investigate the stability of the 
solutions in the different regions of the parameter space, 
we have used the Nyquist criterion. The description of the 
stability diagram is similar to the one given in Sec. 16.3 
and the different possible cases can be understood directly 
from the reading of Figs. [19] and [20] The best way is to 



fix y and progressively increase the value of rj. For y < y„, 
the distribution has only one maximum at Vb = V— so 
the Nyquist curve has one intersection with the x-axis (in 
addition to the limit point (1, 0)). For rj — > 0, we find that 
e r (V_) > so the system is stable. As we increase 77 and 
pass above the solid line, we find that e r (V-) < so the 
system is unstable. For y > y*, the distribution has two 
maxima at Vq = V- and Vq = V+ and one minimum at 
Vq = Vp- At each intersection with a marginal line, one of 
the values e r (V-), e r (V p ) or e r (V+) changes sign. We have 
indicated by symbols like ( — h +) the respective signs of 
e r (V-), e r (V p ) and e r (V + ). We can then easily draw by 
hands the corresponding Nyquist curve. Therefore, it is 
not necessary to show all the possibilities and we have 
only indicated a few representative cases in Figs. [2T1I241 
for illustration. 




Fig. 19. Stability diagram of the asymmetric double- 
humped distribution for A < A* (specifically A = 1.02). 
There exists three marginal branches 77c \y), Vc P \y) & n d 
Vc + \y)- The symbols have been defined in the text. 
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Fig. 20. Stability diagram of the asymmetric double- 
humped distribution for A > A* (specifically A = 10). 
There exists two marginal branches ?7c \y) and r)i + \y). 
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Fig. 21. A < A*: Nyquist curve for y* < y < y c and 
r\ < rji~' (specifically A = 1.02, y = 1.1 and tj = 1). The 
DF has a global maximum at V-, a minimum at V p and 
a local maximum at V+. The DF is stable (with respect 
to this perturbation) because the Nyquist curve does not 
encircle the origin. Case (+ + +). 
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Fig. 22. A < A*: Nyquist curve for y* < y < y c and 
Tji < rj < r]i + ^ (specifically A = 1.02, y = 1.1 and rj = 
3.436). The DF has a global maximum at V_, a minimum 
at V p and a local maximum at V+. The DF is unstable 
(with respect to this perturbation) because the Nyquist 
curve encircles the origin once. There is N = 1 unstable 
mode. Case (- + +). 



Remark: In Figs. Q32 and [501 we see that the system 
is stable for perturbations with k > (fc c )* (correspond- 
ing to the solid line) and unstable for perturbations with 
k < (fc c )*. This critical wavenumber (fc c )* corresponds to 
a marginal perturbation where the phase velocity ui/k co- 
incides with the global maximum vq = V— of the velocity 
distribution. 



8. The case of plasmas 

We shall now compare the previous results obtained for 
self-gravitating systems to the case of plasmas. 
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Fig. 23. A < A*: Nyquist curve for y* < y < y c and 
77c + ' < f] < r]c (specifically A = 1.02, y = 1.1 and 77 = 
4). The DF has a global maximum at V—, a minimum 
at Vp and a local maximum at V+. The DF is unstable 
(with respect to this perturbation) because the Nyquist 
curve encircles the origin twice. There are N = 2 unstable 
modes. Case (- + -). 
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Fig. 24. A < A*: Nyquist curve for y* < y < y c and 
V > vi P) (specifically A = 1.02, y = 1.1 and rj = 6). The 
DF has a global maximum at V- , a minimum at V p and a 
local maximum at V+. The DF is unstable (with respect 
to this perturbation) because the Nyquist curve encircles 
the origin once. There is N = 1 unstable mode. Case (- - 
-)■ 



8.1. A brief historic 

The name plasma was introduced by Tonks & Langmuir 
(1929) to describe an ionised gas made of electrons and 
ions. They found that cold plasmas oscillate with a nat- 
ural pulsation lo p = (4npe 2 /m 2 ) 1 / 2 . In order to take into 
account thermal effects, plasmas were initially described 
in terms of hydrodynamic equations. A dispersion re- 
lation was derived by Thomson & Thomson (1933) for 
isothermal perturbations and by Gross (1951) for three 
dimensional adiabatic perturbations. However, fluid equa- 
tions do not provide a good description of plasmas and 
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rely on arbitrary assumptions concerning the perturba- 
tions (see van Kampen 1957). Indeed, plasmas are essen- 
tially collisionless so that a hydrodynamic description is 
not clearly justified. Vlasov (1938,1945) was the first to 
attempt to derive the dispersion relation directly from 
the collisionless Boltzmann equation (nowdays called the 
Vlasov equation) . He obtained the expression of the pulsa- 
tion u>(k) in the high wavelength limiil 30 !. However, Landau 
(1946) criticized his mathematical treatment showing that 
there are serious divergences in the integrals considered by 
Vlasov. Landau performed a rigorous mathematical study 
using an appropriate contour of integration in the com- 
plex plane and showed that the plasma undergoes damped 
oscillations. The pulsation is in agreement with the ex- 
pression given by Vlasov but, in addition, the plasma 
undergoes collisionless damping (nowdays called Landau 
damping). Following Landau's seminal work, several au- 
thors studied the stability of a plasma. They showed 
that single-humped distributions are always stable (Bcrz 
1956, Appendix by W. Newcomb in Bernstein 1958, Auer 
1958, Penrose 1960, Noerdlinger 1960, Jackson 1960). 
Then, they considered double humped distributions mod- 
eling two contrastreaming beams (Haeff 1949, Nergaard 
1948, Pierce & Hebenstreit 1949) and determined partic- 
ular conditions under which such distributions are unsta- 
ble (Twiss 1952, Buneman 1958, Auer 1958, Kahn 1959, 
Buncman 1959, Penrose 1960, Noerdlinger 1960, Jackson 
1960). They found qualitatively that the plasma becomes 
unstable when the stream velocity becomes sensibly larger 
than the thermal velocity. 

After briefly recalling some classical results, we will 
provide a detailed study of the stability /instability of a 
double-humped distribution made of two Maxwellians. We 
shall complete the study of a symmetric double-humped 
distribution (by explicitly computing the range of unstable 
wavelengths) and consider for the first time the case of 
asymmetric double-humped distributions. We will show 
that the nature of the problem changes above a critical 
asymmetry A* = 3.3105. 

8.2. The Euler-Poisson system 

Wc consider a plasma made of electrons with mass m and 
charge — e and ions with mass rrii and charge +e. We as- 
sume the ions to be infinitely massive so that they do not 
contribute to the motion and only provide a neutralizing 
background. The fluid equations describing the motion of 
the electrons can be written 



dp 
dt 



+ V • {pa) = 0, 



-j£ + (u- V)u = — Vp- V$, 

at p 

Aire 2 

A$ = riP- Pi)- 



(233) 



(234) 



(235) 



The same result was derived in a different manner by Bohm 
& Gross (1949) and Jackson (1960). 



These equations differ from the Euler-Poisson system (JTJ) — 
fl3J) describing a self-gravitating barotropic gas only in the 
sign of the interaction and in the presence of a neutraliz- 
ing background. This neutralizing background avoids the 
Jeans swindle since an infinite homogeneous distribution 
of electrons (p = est, u = 0, $ = 0) is a steady state of the 
fluid equations (|233|) - (|235|) provided that the condition of 
electroncutrality p — pi is satisfied. 

Linearizing the Euler-Poisson system around an infi- 
nite homogeneous distribution of electrons and decom- 
posing the perturbations in normal modes of the form 
e i(k-r-wt)^ we obtain the dispersion relation (Nicholson 
1992): 

lj 2 = c 2 s k 2 + uj%, (236) 
where c 2 = p'(p) is the velocity of sound and 



Airpe 2 



(237) 



is the plasma pulsation. For an isothermal equation of 
state p — pT, we get 



Tk 2 + ui 2 



(238) 



and for a polytropic equation of state p = Kp 1 , we obtain 



iTk 2 + ul 



(239) 



The dispersion relation (|236|) shows that a plasma de- 
scribed by fluid equations is always stable. For k = 0, 
the pulsation tends to a constant 



(240) 



Therefore, in this limit of small wavenumbers, plasma 
oscillations do not propagate in space: p(r, t) = 
p(r,0)e~ lujpt . This behavior is in marked contrast with 
all usual types of oscillating phenomena such as acous- 
tic waves for which uj = ck. For large wavenumbers, 
the dispersion relation approaches that of sound waves: 
u) = ±c s k. 

Remark: Since / SpSQdr = f{A5<f>) 2 dr > for 
a plasma, the second variations (|12[) of the energy func- 
tional (|TU)) arc always positive: 5 2 W > 0. Therefore, the 
energy functional has a unique (global) minimum and the 
plasma is nonlinearly dynamically stable with respect to 
the barotropic Euler-Poisson system. 

8.3. The Vlasov-Poisson system 

The fluid equations are not very appropriate to describe 
a plasma because most plasmas are in a regime where 
the collisions between charges are negligible. Therefore, 
from now on, we shall describe Coulombian plasmas by 
the Vlasov-Poisson system 



df df df 
— + v • — — V$ • — 
dt dr dv 

47re 2 

A$ = rip- Pi), 



0. 



(241) 



(242) 
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which ignores collisions between charges. These equations 
are similar to the Vlasov-Poisson system (|3"3")) - (f3~4")) describ- 
ing collisionless stellar systems except that the sign of 
the interaction is reversed and that there is a neutraliz- 
ing background. We shall investigate the linear dynamical 
stability of an infinite homogeneous medium. Therefore, 
we consider stationary solutions of the Vlasov equation of 
the form / = /(v). The dispersion relation can be written 
(Nicholson 1992): 



e(k,cj) 



1 



Ane 2 



m- 



/'(«) 



dv = 0. 



(243) 



c 



As in Sec. 13.31 we have taken k along the z-axis and noted 
v for v z and f(v) for J f(\)dv x dv y . 

In general, the dispersion relation (|243[) cannot be 
solved explicitly to obtain oj(k) except in some very sim- 
ple cases. For example, for cold systems described by the 
distribution function / = pS(v — vq), we obtain after an 
integration by parts 



LO = V()k ± Ul p . 

In particular, when vq = 0, we get 



(244) 



(245) 



The system is stable to all wavenumbers. We also note that 
the dispersion relation (|245|) coincides with the dispersion 
relation (|236|) obtained with fluid equations with c s = 0. 

Let us briefly recall important properties of the disper- 
sion relation (|243[) obtained in plasma physics (Balescu 
1963, Nicholson 1992). We consider a symmetric distri- 
bution and write the complex pulsation in the form ui = 
uj r + iuji. When u>i <C w r (weakly damped perturbations), 
the complex pulsation is given bjJ3 



ZTk 2 



2pk 2 \ k 



(246) 



with T = (v 2 ) (where we recall that v = v z in the present 
case). These relations are valid for fc/fco -C 1. The real 
part of the pulsation satisfies a dispersion relation of the 
form (|B.17j) with c 2 = 3T. Therefore, large wavelengths 
perturbations in a collisionless plasma correspond to one 
dimensional iscntropic perturbations with index 7 = 3 in 
a gas (see Appendix |B|) . This relation dispersion is called 
the Langmuir wave dispersion relation (Nicholson 1992). 
On the other hand, when f'(uj p /k) < 0, the expression of 
the imaginary part of the pulsation corresponds to Landau 
damping. If f'(v) < for all v > 0, as for the Maxwellian, 
the plasma is stable. For the Maxwell distribution, we get 
the Landau formula 



kp 
k 



(247) 



31 A more accurate expression (Jackson 1960, Nicholson 1992) 
is obtained by replacing f'(uj p /k) by f'(uj r /k). This introduces 
a factor e 3 ^ 2 in the Landau formula (124711. 



for the damping rate, where we have introduced the Debye 
wavenumber 



a _ 47re 2 ffp 2 
k D = — = /3uj 



(248) 



where (3 = 1/T is the inverse temperature. Alternatively, 
if there exists wavenumbers k such that /'(w p /fc) > 0, 
the perturbation grows and the plasma is unstable. Since 
the expression (|246[) is valid for small k, this implies that 
f'(v) must be positive for large velocities. This corre- 
sponds for example to the "bump-on-tail" situation ana- 
lyzed in plasma physics (Nicholson 1992). The dispersion 
relation ([Mo]) -a was first obtained by Vlasov (1938,1945). 
However, the derivation given by Vlasov was not rigorous 
since he evaluated the integral (|243|) on the real axis on 
which there is a singularity. Landau (1946) performed a 
rigorous mathematical study and showed that, in addition 
to the oscillations, the waves must be damped exponen- 
tially leading to equation (|246p ~b for the damping rate. It 
is interesting to note that the Langmuir wave dispersion 



relation (|246[) -a is rather independent on the precise form 
of the distribution function (it depends only on the vari- 
ance T — (v 2 )) while the Landau damping (|246[) -b is very 
sensitive to the distribution function. For example, the ex- 
pression of Landau damping for polytropic distributions 
is given by (see Appendix E of Chavanis & Dclfini 2009): 



1 - 



1 



1 



n , 

2 1 n+l 



k 2 



2(n + l)fc 2 



ra-3/2 



(249) 



When u r <C (heavily damped perturbations), the 
pulsation is given for the Maxwellian distribution by 

7TWr, k 2Ll> v , r - 



±971- 



(250) 



These expressions are valid for k/kjj +00. There ex- 
ists several branches of solutions parameterized by the 
odd integer m (but the branch m = 1 is the most rele- 
vant). Expression (|250|) of the complex pulsation was first 
obtained by Landau (1946). The evolution of u> r (k) and 
LOi{k) with k for a Maxwellian distribution are represented 
in Fig. 3 of Jackson (1960) and it can be compared with 
Fig. [21 for stellar systems. For small k, the damping is 
weak but it increases rapidly for increasing wavenumbers 
so that the waves are never observed at sufficiently large 
wavenumbers. 

Remark: Let us consider a plasma with a DF of the 
form / = f(v 2 ) with /' < 0. Since J SpSQdr > 
for a plasma, the second variations (|3"7]l of the con- 
strained pscudo entropy (|3"5T) are always negative: S 2 S < 0. 
Therefore, the pseudo entropy has a unique (global) max- 
imum at fixed mass and energy and the plasma is nonlin- 
early dynamically stable. These variational technics were 
introduced in the plasma literature by W. Newcomb in 
Bernstein (1958) and Gardner (1963); see also the review 
of Holm et al. (1985). 
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8.4. The condition of marginal stability 

When uji = 0, the real and imaginary parts of the dielectric 
function e(k,u> r ) = e r (k,uj r ) + iei{k,u> r ) are 



e r (k,uj r ) = 1 



ei(k,u) r ) 



m 2 k 2 
An 2 e 2 



m 



■k 2 



/'(Wr/ft). 



(251) 



(252) 



The condition of marginal stability corresponds to 
e(fc,w) = and Wj = 0, i.e. e r (fc,w r ) = ei(k,co r ) = 0. 
We obtain therefore the equations 



m 2 k 2 



v — uj r /k 



(253) 



Furthermore, repeating the steps of Sec. 
k -> k r thafi 



47r 2 e 2 /"(0) 



1 fc 2 



(fc -> fc c ). 



we obtain for 



(260) 



This formula leads to the following result. If the distribu- 
tion is maximum at v = 0, so that /"(0) < 0, we find that 
the mode u = iu>i is stable for k < k c and unstable for 
k > k c . Alternatively, if the distribution is minimum at 
v = 0, so that /"(0) > 0, we find that the mode u — itOi is 
stable for k > k c and unstable for k < k c . This result will 
be illustrated in connection to Fig. [^Hl for the symmetric 
double humped distribution. Note that this result implies 
that a symmetric distribution satisfying inequality (|258[) 
is always unstable (to some wavenumbcrs) . 



/' (Ur/k) = 0. 



(254) 



The second condition (|254l) imposes that the phase ve- 
locity uj r /k = v ex t is equal to a velocity where f(v) is 
extremum (f'(v ex t) = 0). The first condition (|253[) then 
determines the wavenumber(s) corresponding to marginal 
stability. It can be written 



47re 2 



f'(v) 

' ~ V ext 



dv 



1/2 



(255) 



Finally, the pulsation(s) corresponding to marginal sta- 
bility are u r = v ex tk c and the perturbation behaves like 
Sf ~ e ~ iUrt . 

8.5. Particular solutions of e(k,u) = 

We can look for a solution of the dispersion relation 
e(fc,w) = in the form lo = iuji corresponding to uj r = 0. 
In that case, the perturbation grows (wi > 0) or decays 
(u>i < 0) without oscillating. Like in Sec. 13. 81 we shall as- 
sume that f(v) is an even function. For cjj > 0, the growth 
rate is given by 



i-i^r^u= . 



m 2 k 2 



(256) 



00 « + # 



For Ui < 0, the decay rate is given by 



47re 2 
m 2 k 2 



+ CJO 



vf'(v) 



dv 



in 2 e 2 



— oo 



1)1 



'■k 2 



lUJj 



f ( _1 ) = o. (257) 



Let us assume that the distribution satisfies 

, + oc f{v) 



■dv>0. 



V 



(258) 



This is not the generic case. For example, the Maxwell 
distribution does not satisfy this inequality. If inequality 
(|258[) is satisfied, then the marginal mode w = corre- 
sponds to the critical wavenumbcr 



47re^ 



dv 



1/2 



(259) 



8.6. The Nyquist method 

To apply the Nyquist methods, we have to plot the curve 
(e r (fc, 0J r ), €i(k, u) r )) parameterized by u> r going from — oo 
to +oo (for a fixed wavenumber k). Let us consider the 
asymptotic behavior for uj r — > ±oo. Since f(v) is posi- 
tive and tends to zero for v —> ±oo, we conclude that 
ti{k,to r ) — > for uj r — > ±oo and that 6i(k,u r ) < for 
uj r — > — oo while ei(k,cu r ) > for cu r — > +oo. On the other 
hand, for w r — > ±oo, we obtain at leading order 



e r (k,uj r ) ~ 1 



47r 2 e 2 p 



(cu r — > ±oo). 



(261) 



From these results, we conclude that the behavior of the 
curve close to the point (1,0) is like the one represented in 
Fig. 1251 For u r /k = v ext , where v ext is a velocity at which 
the distribution is extremum (f'{v ex t) = 0), the imaginary 
part of the dielectric function ej(fc, kv ex t) = and the real 
part of the dielectric function 



/, , s 47r 2 e 2 f 00 f'iv) 
e r (k, kv ext ) = 1 — / dv. 



i 2 k 2 



(262) 



Subtracting the value f'(v ext ) = in the numerator of the 
integrand, and integrating by parts, we obtain 

e r (k,kv ext ) = H 57^/ — ^— dv. (263) 



m 2 k 2 



{v - v ext y 



If VMax denotes the velocity corresponding to the global 
maximum of the distribution, we clearly have 

e r {k 1 kv M ax)> I- (264) 

32 In the derivation, we have assumed that f"(0) =fc 0. If 
/"(0) = 0, we need to develop the Taylor expansion to the 
next order. 

33 Nyquist (1932) introduced his graphical method in rela- 
tion to servomechanism theory. His method was first applied to 
plasma physics by Harris (1959), Penrose (1960) and Jackson 
(1960). 
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8.6.1. Single-humped distributions 

Let us assume that the distribution f(v) has a single max- 
imum at v = w (so that /'(i>o) = 0) and tends to zero 
for v — > ±oo. Then, the Nyquist curve cuts the x-axis 
(ei(k,uj r ) vanishes) at the limit point (1,0) when ui r — > 
±oo and at the point (e r (k, kvo), 0) when uj r /k = Vq. Due 
to its behavior close to the limit point (1, 0), the fact that 
it rotates in the counterclockwise sense, and the property 
that e r (k,kvo) > 1 according to Eq. (|264[) . the Nyquist 
curve must necessarily behave like in Fig. 1251 Therefore, 
the Nyquist curve starts on the real axis at e r {k,uj r ) — 1 
for ui r — > — oo, then going in counterclockwise sense it 
crosses the real axis at the point e r (k,kvo) > 1 and re- 
turns on the real axis at e r (k,ui r ) = 1 for u> r — > +oo. 
Therefore, it cannot encircle the origin. According to the 
Nyquist criterion exposed in Sec. 13.51 we conclude that a 
single-humped distribution is always linearly dynamically 
stableR 
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Fig. 25. Nyquist curve for the Maxwell distribution. The 
DF is always stable. 



Let us specifically consider the Maxwell distribution 
(|102p for illustration. The dielectric function (|243[) associ- 
ated to the Maxwell distribution is 



Introducing the Debye wavenumber (|248[) . it can be 
rewritten 



. , kn ( \fBu) 

e(k,u) = 1 + yf W' 



k 2 \ k 



(266) 



It will be convenient in the following to work with di- 
mensionless quantities. We introduce the dimensionless 
wavenumber and the dimensionless pulsation 



(267) 



The fact that a single-humped distribution is always stable 
is known as Gardner (1963)'s theorem; see also Jackson (I960). 



Noting that y/]3u)/k = y^O, the dielectric function (|266l) 
can be rewritten 



iri,n) = l + r)W(SfjSl). 



(268) 



When Qi = 0, the real and imaginary parts of the di- 
electric function e(rj, fl r ) = e r (?7, ^r) + itiij], ^r) are given 

by 



e r (r), n r ) = 1 + T)W r (y/rjtt r ) ■ 
£i(r/, fJ r ) = r\Wi (y/rjttr) ■ 



(269) 
(270) 



The Nyquist curve for the Maxwell distribution is repre- 
sented in Fig. l25l 

On the other hand, if we consider a polytropc of index 
n = 1/2 (waterbag distribution) as in Sec. 15.51 we find 
that the dielectric function is given by 



kn I OJ 



kVT 



(271) 



where Wi/2(x) is given by Eq. (|175[) . The condition 
e(fc, ui) = determines the dispersion relation. We get 



lo 2 =luI + 3Tk 2 . 



(272) 



We note that, for the water-bag distribution, the general 
asymptotic behavior (|246[) becomes exact. We also note 
that for the specific index n = 1/2 (7 = 3), the dispersion 
relation in a collisionless stellar system takes the same 
form as in a gas (see Sec. 18. 2p . 

8.6.2. Double-humped distributions 

Let us consider a double-humped distribution with a 
global maximum at VMax, a minimum at v m i n and a local 
maximum at v max . We assume VMax < «min < v max . The 
Nyquist curves starts at (1,0), progresses in the counter- 
clockwise sense and crosses the rr-axis at e r (vMax) > 1, 
then at e r (v m i n ) and e r (v max ). We can convince ourselves 
by making drawings of the following results. If 

(+ + +): e r (v Ma x) > 0, e r (v mm ) > 0, e r (v max ) > 0, 

(H ): e r (v M ax) > 0, e r (v mm ) < 0, e r (v max ) < 0, 

(+ + -): £ r (vMax) > 0, e r (v mm ) > 0, e r (v max ) < 0, 
the Nyquist curve does not encircle the origin so the sys- 
tem is stable. If 

(H h): £ r (vMax) > 0, e r (v min ) < 0, e r ( )>o, 

the Nyquist curve rotates one time around the origin so 
that there is one mode of instability. Since e r (vMax) > 
there is no mode of marginal stability with uj r /k = VMax- 

Cases (+ + +), (H h) and (H ) are observed in Sec. 

18.81 for an asymmetric doublc-humpcd distribution made 
of two Maxwellians. 

If the doublc-humpcd distribution is symmetric with 
respect to the origin with two maxima at ±i>* and a min- 
imum at v = 0, we get the same results as above with the 
additional properties e r (vMax) = £r{vmax) = £r(v*) > 1 
and e r (v m i n ) = e r (0). We have only two cases (+ + +) 
and (H h). They are observed in Sec. 18.71 for a symmet- 
ric double-humped distribution made of two Maxwellians. 
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Since e r («*) > 0, there is no mode of marginal stability 
with LJ r /k = ±u*. 

It can be shown that a plasma is unstable (to some 
wavelengths) iff f(v) has a minimum at a value v = v m i n 
such that 



f'(v) 



-dv 



v - v„ 



f(v) - f{v m in) 
[V — Vmin ) 



dv > 0. (273) 



This result was proven by Penrose (1960) and it is some- 
times called the Penrose criterion^!. This criterion can be 
deduced from the Nyquist method as follows. A double- 
humped distribution f(v) is unstable if there exists a 

range of k such that we are in the situation (H h), i.e. 

£r(v m in) < and e r (v max ) > 0. The first condition can 
be satisfied (for sufficiently small k) iff condition (|273D is 
realized, which is the Penrose criterion. Then, the range 
of unstable wavenumbers is 



(274) 

If the first integral is negative, the range of unstable 
wavenumbers is k < kc ■ This corresponds to Fig. 
|3"T1 and to the region y > y c in Fig. [3U If the first in- 
tegral is positive, the range of unstable wavenumbers is 
^(jnax) ^ ^ ^ This corresponds to the region 

y* < y < y c in Fig. [32] (indeed, y c corresponds precisely to 
the case where the first integral becomes equal to zero). 



8.7. The symmetric double-humped distribution 

Let us now consider the symmetric double-humped distri- 
bution (|184p . In the plasma case, the dielectric function 
can be written 

e( v , n) = 1 + 1 [W(^rjn - Vy)) + w(^n + ^y))] , 

(275) 

where rj (wavelength) and (pulsation) are defined by 
Eq. (f2"6T]) and y (separation) is defined by Eq. (fTHSj) . The 
condition of marginal stability corresponds to e{k,u) = 
and uji = 0. The condition €i(k,u) r ) — is equivalent 
to f'(ui r /k) — so that the phase velocity ui r /k = vq 
corresponds to the velocities where the distribution is ex- 
tremum. Then, the wavenumbers at which the distribution 
is marginally stable are obtained by solving e r (k,u r = 
kvo) = 0. Proceeding as in Sec. [S] and introducing the 
parameter x = Vo/V a , we find that the equations deter- 
mining the critical wavenumbers rj c (y) are given by 



y 



i 

2x 



In 



1 



1 - x 



[W r (^y(x-1)) + W r (^y(x + 1))] 



(276) 



(277) 



We note that only the sign in Eq. (|277p changes with re- 
spect to the study of the gravitational case, so we can 
readily adapt the results of Sec. [S] to the present situation 
by simply reverting the sign. For ^JrjSlr = ±T4, corre- 
sponding to x 7^ 0, there is no physical solution to Eqs. 
(I276l) - (|277|) with positive rj > 0. Therefore, in the plasma 
case, there is no marginal mode with non zero pulsation for 
the symmetric double-humped distribution (in agreement 
with the general discussion of Sec. 18.6.2"]) . We now consider 
the marginal mode with O r = 0. This corresponds to the 
"degenerate" solution x = (for any y) for which Eqs. 
(|276 > (|2"77 ]) reduce to 



-l 



(278) 



According to Fig. physical solutions exist only for y > 
Vmax = = 1.708. We note that the range of parameters 
that was forbidden in the gravitational case is now allowed 
in the plasma case and vice versa. For y < y ma x, the 
system is stable. This result is to be expected since, for 
y = 0, the distribution (|184[) reduces to the Maxwellian 
that is stable for a repulsive interaction. We now consider 
the range y G [y ma x, +oo[. For y -> y max , we have 



vi 0) (y) 



]J Umax 



(y ymax)- 



On the other hand, the curve (y) has a minimum 
at (4.511,3.512) (see Appendix A of Chavanis & Delfini 
2009). Finally, for y — > +00, we have 

(0) (y) 



(279) 



y 



+00, 



(y 



(280) 



Again this result is expected because, for y — > +00, the 
two humps do not "see" each other and behave as two 
independent single-humps distributions that are stable in 
the repulsive case. 

The critical wavenumber rji ^ (y) corresponding to 
marginal stability determined previously is represented as 
a function of the separation y in Fig. [25] We have also 
plotted the line y = 1. On the left of this line, the dis- 
tribution has a single maximum at Vq = and on the 
right, the distribution has two maxima at Vq = ±V* and 
a minimum at Vq = 0. In order to investigate the stability 
of the solutions in the different regions, we have used the 
Nyquist criterion. 

For y < 1, the DF has only one maximum at Vo = 0. 
There is no marginal mode and the distribution is always 
stable (see Fig. [27]). 

For 1 < y < y m ax, the DF has a minimum at Vo = 
and two maxima at ±V*. There is no marginal mode and 
the distribution is always stable (see Figs. |2"51) . 

For y > y m ax, the distribution has a minimum at Vo = 
and two maxima at ±14. There exists one wavenumber 
r)c 0%> at which the distribution is marginally stable. For 
rj = rjc , the marginal perturbation does not propagate 
(f2 r = 0). By considering the Nyquist curves in this region 



In fact, we found that an equivalent criterion was obtained ( see Figs. 



we find that the DF is stable for 77 < rj, 



(0) 



at the same period by Noerdlinger (1960). 



and unstable for rj > rj. 



(0) 
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Fig. 26. Stability diagram of the symmetric double- Fig . 2 8. Nyquist curve for 1 < y < y max (specifically 
humped distribution (HMJ) in the case of plasmas. y = 12 an d 77 = 6). The DF is stable for any perturbation. 

Case (+ + +). 
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Fig. 27. Nyquist curve for y < 1 (specifically y = 0.5 and 
77 = 6). The DF is stable for any perturbation. Case (+). 



Let us make some remarks: 

1. If we consider all possible perturbations, Fig. [25] 
shows that a symmetric double-humped distribution is 
stable if y < y ma x (i-e. v a < 1.307 T 1 / 2 ) and unstable 
if V > Umax (i-e. v a > 1.307 T 1/2 ). Therefore, instabil- 
ity occurs when the drift velocity v a is sensibly larger 
than the thermal speed y/T so that the two humps are 
well separated. In particular, at T = 0, a double-humped 
distribution is always unstable to some wavelengths (see 
Appcndix|E|. More precisely, a symmetric double- humped 
distribution with y > y max is stable for perturbations 
with wavenumbers k > (fc c )* (corresponding to the solid 
line) and unstable for perturbations with wavenumbers 
k < (fc c )*. The critical wavenumber (k c )* corresponds to 
the marginal perturbation for which the phase velocity 
uj/k coincides with the minimum of the velocity distribu- 
tion: Vq = 0. 

2. For y > y maX) the mode u> = iuji is stable for r\ < r/i ^ 
and it becomes unstable when we increase 77 above the 
marginal line 77 = i]^ . This is consistent with the general 
result (|260[) since, for y > 1, the DF is minimum at vq = 0. 



oT 0- 




Fig. 29. Nyquist curve for y > y max and 77 < r]^ (specif- 
ically y — 5 and 77 = 2). The DF is stable for this range of 
wavenumbers. Case (H — | — h). 



3. For a double- humped distribution, the critical wave- 
length is infinite for v a < \Jy m axT (stable), then de- 
creases, reaches a minimum and increases again. By con- 
trast, if the plasma is modeled as a contrastrcaming gas, 
the critical wavelength is infinite for v a < c s (stable), then 
jumps discontinuously to zero at v a = c s and increases for 
v a > c s (Ikeuchi et al. 1974). 



8.8. The asymmetric double-humped distribution 

The dielectric function associated to the asymmetric 
double-humped distribution (|207[) in the plasma case is 



efaSl) = 1 + 



n 



W{^jn - y/y)) 



1 + A 

+Aw(^jn + x /j])) 



(281) 
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Fig. 30. Nyquist curve for y > y max and r) > r]^ (specif- 
ically y = 5 and rj = 6). There is one mode of instability 
in this range of wavenumbers. Case (- + -). 







< x c) 


A < A, 

(+ - +) 






\ «,= V P 


Unstable ^ 
(one mode) 


- 

(+) 












(+ + +) 


Stable 


y. 









Fig. 31. Stability diagram of the asymmetric double- 
humped distribution for A < A* (specifically A = 1.02). 
There exists one marginal branch ry p (y) ■ 



Proceeding as in Sec. [7] and introducing the parameters 
x = V~o/V a and y = V a 2 , we find that the equations deter- 
mining the critical wavenumbers rj c (y) are given by 



y 



v 



2x 



1 



ln(A) 
2x 

1 + A 



[W r (^y(x - 1)) + AW r (^y(x + 1))] 



(282) 



(283) 



Equation (|282p determines the extrema of the distribution 
f(v) and Eq. 



(|283|) determines the wavenumbers corre- 
sponding to the modes of marginal stability. As in Sec. 
the curve rj c (i/) can be obtained by varying x between —1 
and +1. In the plasma case, there exists physical solutions 
with positive rj only for < x < x c . For x — > + , using 
, we find that y — > +oo. Then, we get 

y->+oo, (y^+oo). (284) 

we find that y — > y c and 

+oo, (y -> y c ), (285) 



Eq. (1282 

For x — > x 



y - Vc 

where s = p if A < A* and s = + if A > A» . For A < A* , 
we get only one marginal branch rji (y) corresponding to 
the mode fl r = V p /^/rfc > (see Fig.[3T]). For A > A*, we 

get two marginal branches r]i P \y) and rjc + \y) correspond- 
ing to the modes £l r — V p / y/T% > and O r — V+/ ^frfc > 
(see Fig. l32j) . They connect each other at corre- 
sponding to x = x*. At that point they touch the line 
y = y„ separating distributions with one or two maxima. 

The stability diagrams corresponding to the asymmet- 
ric double-humped distribution with A < A* and A > A* 
are represented in Figs. [31] and [32j To investigate the 
stability of the solutions in the different regions of the 
parameter space, we have used the Nyquist criterion. If 
A < A*, the description of the stability diagram can be 
made in continuity with the case of a symmetric double- 
humped distribution: for y < y c (i.e. v\ < y c (A)T), the 
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Fig. 32. Stability diagram of the asymmetric double- 
humped distribution for A > A* (specifically A = 50). 
There exists two marginal branches r]c P \y) and r]c + \y). 



plasma is stable to all perturbations; for y > y c (i.e. 
v a > 2/c(A)T), the plasma is stable for k > k*f\v a ) 
and unstable for k < kc P \v a ). If A > A*, the stabil- 
ity diagram changes: for y < y* (i.e. v 2 < j/*(A)T), the 
plasma is stable to all perturbations; for y* < y < y c 
(i.e. y*(A)T < v\ < y c (A)T), the plasma is stable for 

k > kc P \v a ), unstable for ki (v a ) < k < kc P \v a ) and sta- 
ble for k < kc (v a ): this is similar to a re-entrant phase; 
for y > y c (i.e. v 2 a > y c (A)T), the plasma is stable for 
k > kc(v a ) and unstable for k < ki P \v a ). In conclusion, 
instability occurs when the drift velocity is sensibly larger 
than the thermal speed so that the two humps are well 
separated. However, the precise threshold changes whether 
the asymmetry is smaller or larger than A*. For A < A*, 
the plasma is stable iff v\ < y c (A)T and for A > A*, the 
plasma is stable iff v 2 < y»(A)T. 
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Fig. 33. A > A*: Nyquist curve for y* < y < y c and 
rj < r]c (specifically A = 50, y = 5 and r\ — 5). The 
DF has a global maximum at V-, a minimum at V p and 
a local maximum at V+. The DF is stable (with respect 
to this perturbation) because the Nyquist curve does not 
encircle the origin. Case (+ + +). 



30 1 — p 




-10 10 20 30 

e 

r 

Fig. 35. A > A*: Nyquist curve for y* < y < y c and 
rj > 77c + ' (specifically A = 50, y = 5 and ?y = 30). The 
DF has a global maximum at VL, a minimum at V p and 
a local maximum at V+. The DF is stable (with respect 
to this perturbation) because the Nyquist curve does not 
encircle the origin. Case (H ). 
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Fig. 34. A > A*: Nyquist curve for y* < y < y c and 
r)c < rj < rji + ^ (specifically A = 50, y = 5 and 77 = 15). 
The DF has a global maximum at V-, a minimum at V p 
and a local maximum at V+. The DF is unstable (with 
respect to this perturbation) because the Nyquist curve 
encircles the origin once. There is N = 1 unstable mode. 
Case (H h). 

9. Conclusion 

We have carried out an exhaustive study of the linear 
dynamical stability/instability of an infinite and homo- 
geneous self-gravitating medium with respect to pertur- 
bations with wavenumber k. This is the classical Jeans 
problem that describes the growth of perturbations and 
the emergence of structures (like stars and galaxies) in 
astrophysics and cosmology. We have considered the case 
of a collision-dominated gas described by the Euler equa- 
tion or the case of a collisionless stellar system described 
by the Vlasov equation. In the latter case, we have stud- 
ied single-humped distributions as well as symmetric and 
asymmetric double-humped distributions. We have used 



the Nyquist theorem to settle the stability of the system. 
Detailed stability diagrams in the (v a ,k) plane for dif- 
ferent values of the asymmetry parameter A have been 
obtained. We have also derived general analytical results 
concerning the dispersion relation. Finally, we have stud- 
ied the case of Coulombian plasmas in parallel. The gen- 
eral methods developed here to study the stability of a ho- 
mogeneous self-gravitating medium can have applications 
in other domains, not necessarily astrophysics, where the 
system is described by the Euler or the Vlasov equations 
coupled to a mean field equation for the potential. One 
example is the HMF model (Chavanis & Delfini 2009). 
Although the structure of the problem is the same, the 
results and the stability diagrams differ because the po- 
tential of interaction is different. Many generalizations of 
our work are possible, changing the potential of interac- 
tion or the class of distribution functions, but the present 
paper provides many methods and analytical results that 
can be useful to investigate more general situations. 

Appendix A: Equivalence between the variational 
problems (gB) and (fill) 

The equivalence between the variational problems (J3SJ) 
and has been discussed in previous works (Chavanis 
2006a, 2008b). We shall here provide a new derivation of 
this result. To that purpose, we shall show the equivalence 
between the stability conditions (fT2j) and ((39)) . 

First of all, using the identities /'(e) = — /3/C"(/) (see 
Sec. [HIE]) and p'(p) = -p/p'($) (see Sec. l2TT]t . we can 
rewrite the second variations (|3"9"|) and (fT2")l in the form 

5 2 F = - l -J { j^dvdw+ l -j5 P 5^dv, (A.l) 

62w= -\S%) dr +\j 5pmdv - (a - 2) 
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We shall determine the perturbation 5/*(r, v) that min- 
imizes S 2 F[Sf] given by Eq. (|A.1|) with the constraint 
Sp = f Sf dv. Since the specification of Sp determines (5$, 
hence the second integral in Eq. (jA.lj) . we can write the 
variational problem in the form 



6{ -2 



(Sf) 2 
f'(e) 



drdv) - / A(r)<5 / Sf dv ) dr = 



(A.3) 



where A(r) is a Lagrange multiplier. This gives 

Sf. = -A(r)/'(e), (A.4) 

and it is a global minimum of S 2 F[Sf] with the previ- 
ous constraint since S 2 (S 2 F) = — | J drdv > (the 
constraint is linear in Sf so its second variations vanish). 
The Lagrange multiplier is determined from the constraint 
Sp = f Sf dv yielding Sp = —A / f'(e)dv. Therefore, the 
optimal perturbation (|A.4[) can be finally written 



Sf, 



Sp 



ff'(e)dv 



f'(e). 



(A.5) 



Therefore, S 2 F[Sf] > S 2 F[5f,}. Explicating 5 2 F[8f*\ using 



Eqs. (jA.lj) and (|A.5|) . we obtain 
S 2 F[Sf] >-\j 



(M 2 dr+ l 



SpSQ dr. 



(A.6) 



Finally, using //($) = f /'(e) dv, the foregoing inequality 
can be rewritten 



2F[5f ^-\l%) dV+ \j Sp5 * dV 



S 2 W[Sp], 



(A.7) 



where the r.h.s. is precisely the functional (jA.2j) . 
Furthermore, there is equality in Eq. (|A.7j) iff Sf = Sf,. 
This proves that the stability criteria (jT2"j) and (j3"9"j) are 
equivalent. Indeed: (i) if inequality (j!2D is fulfilled for all 
perturbations Sp that conserve mass, then according to 
Eq. (jA.7j) . we know that inequality (f3T)j) is fulfilled for all 
perturbations Sf that conserve mass, (ii) if there exists a 
perturbation Sp* that makes (5 2 >V[<5p] < 0, then the per- 
turbation Sf, given by Eq. (| A.5|) with Sp = Sp, makes 
S 2 F[Sf] = S 2 W[Sp] < 0. In conclusion, the stability crite- 
ria (fT2"j) and (j3T)j) are equivalent. 

Appendix B: Fluid equations 

In a collisional gas, the evolution of the distribution func- 
tion /(r, v, t) is governed by Boltzmann's transport equa- 
tion 



d -l + ,. d -l-^. d -l^( d -l 
dt dr dv V dt 



(B.l) 



coll 



where (df/dt) co u is the collision operator. This opera- 
tor locally conserves the mass, the impulse and the ki- 
netic energy. The Boltzmann equation relaxes towards the 
Maxwellian distribution for t — > +oo due to collisions. 



From Eq. (jB.lj) , we can easily derive a hierarchy of hydro- 
dynamic equations satisfied by the local moments of the 
velocity distribution. Defining the density p = f f dv, the 
local velocity u = j fv dv and the kinetic temperature 
| pT = | f fw 2 dv (where w = v — u is the relative veloc- 
ity and d the dimension of space) , the first three equations 
of this hierarchy are 



dp 
di 

duj 

at 



V • (pu) = 0, 



dm 



1 OF, 



dxj p dxj 



9$ 



dT 
~dt 



u- VTJ 



Vq 



Pi — 

11 dxi 



(B.2) 



(B.3) 



(B.4) 



where Pij = f fwiWj dv is the pressure tensor and q = 
^ f fw 2 w dv is the current of heat. These equations are 
known as the Maxwell equations of transfer (van Kampcn 
1957) or as the Jeans equations (Binney & Trcmainc 
1987). For a collision-dominated gas, we can close this hi- 
erarchy of hydrodynamic equations by using a Chapman- 
Enskog expansion which is valid when the mean free 
path in the gas is short. The zeroth-order approximation 
amounts to making the local thermodynamic equilibrium 
(L.T.E.) assumption 



fL.T.E.(r,v,t) 



1 



2irT(r,t) 



■ (v-u(r,t)r 

p(r,t)e 2T ^^ 



(B.5) 



This yields q = and = pSij where P = h J fw 2 dv = 
pT is the pressure. In that case, we obtain the Euler equa- 
tions (Huang 1966) 



^+V-(pu) = 0, 

— + u • Vu = — Vp - V$. 

dt p 



dT 
~dt 



u • VTJ + TV • u = 0. 



P = pT- 

Equation (jB.8j) can be rewritten in the form 



d_ 

dt 



P 



T d/2 



0, 



(B.6) 

(B.7) 

(B.8) 
(B.9) 

(B.10) 



where d/dt = d/dt + u • V is the material derivative. This 
equation expresses the conservation of the entropy by the 
Euler equation (dissipativc effects appear at the next order 
in the Chapman-Enskog expansion leading to the Navier- 
Stokcs equations). We have therefore a closed system of 
hydrodynamic equations. The steady states correspond to 
u = and to the condition of hydrostatic equilibrium 
©. However, at the level of the Euler equations, nothing 
determines T(r) at equilibrium, i.e. there exists steady 
states of Eqs. (|B.6j) - (|B.9p with any temperature profile 
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T(r). Two cases are particularly relevant. If we consider 
that the collisions had time to establish a statistical equi- 
librium state, then T(r) = T is uniform at equilibrium 
yielding an isothermal equation of state p(r) = p{v)T. 
By contrast, if collisions had not time to establish a sta- 
tistical equilibrium state, we can consider a steady state 
where the specific entropy is uniform in the whole sys- 
tem: s(r) = s. This is equivalent to having p(r)/T(r) d / 2 , 
p(r)/T(r)( d + 2 )/ 2 or p(r)/ p(r)^ d+2 ^ d spatially uniform at 
equilibrium. We then obtain a polytropic equation of state 
p(r) = Kp(r) 1 with 



d+2 



7 



(B.H) 



In that case, the temperature behaves like T(r) = 
Kp(r)' y ~ 1 . The isothermal distribution is justified by tak- 
ing into account dissipative effects (e.g. using the Navier- 
Stokes equations) and the polytropic distribution is jus- 
tified by completely ignoring dissipative effects. However, 
even if an isothermal distribution has been established un- 
der the effect of collisions, when we investigate the dynam- 
ical stability of such a distribution, we can usually ignore 



dissipative effects and use the Euler equations (|B.6|) - (|B.9p . 
This is because the dynamical time is in general much 
shorter than the collisional relaxation time. 

If we consider the stability of an infinite homogeneous 
self-gravitating system or plasma described by the hydro- 
dynamic equations (|B.6[) - (|B.9|) . we find that the dispersion 
relation is of the form (fT3")) or (|236[) with 



dp 



(B.12) 



In the present case, Sp/Sp ^ p'(p) since the gas is 
not barotropic. Linearizing Eq. (|B.I0|) around the steady 
state, it reduces to 



d 



(B.I3) 



so that p/T d / 2 , p /T^ +2 ^ 2 or p/p( d + 2 )/ d arc independent 
on time in the linear regime. Therefore, the perturbations 
arc adiabatic and wc have 



dp _ d + 2Sp 
p dp' 

yielding Sp = -fTSp hence 
c 2 = 7 T. 



(B.14) 



(B.I5) 



Regrouping all these results, the dispersion relation for a 
self-gravitating system or a plasma described by the (non 



barotropic) Eulcr equations (|B.6|) - (|B.9[) arc 

u? = -/Tk 2 - A-KGp, 



(B.17) 



The are valid for any equilibrium distribution, e.g. isother- 
mal or polytropic. 



Let us now consider the case of a collisionless system 
described by the Vlasov equation obtained from Eq. (|B.I[) 
by neglecting the collision term. The hydrodynamic equa- 
tions (|B.2j) - (|B.4j) remain valid but, in the present case, the 
L.T.E assumption is not justified anymore so that the hier- 
archy of equations cannot be closed. However, for k — > 0, 
the dispersion relation based on the Vlasov equation is 
given for a self-gravitating system by (see Sec. 



AirGp - 3Tfc 2 



(B.I8) 



This dispersion relation is consistent with Eq. (|B.16|) with 
7 = 3 (corresponding to d — I). Therefore, for large wave- 
lengths, a collisionless stellar system can be treated by 
fluid equations experiencing a one dimensional adiabatic 
perturbation. For k — > 0, the dispersion relation based on 
the Vlasov equation are given for a plasma by (see Sec. 



lu 2 = w 2 + 3Tfc 2 



p ft 



r 



2pk 2 J V k 



(B.19) 



The Landau damping is due to collective effects and 
cannot be obtained in a hydrodynamic approximation. 
However, the dispersion relation for the pulsation to r is 
consistent with Eq. (IB. 1 7|) with 7 = 3 (corresponding to 
d = I). Therefore, for large wavelengths, the oscillations 
of a collisionless plasma can be treated by fluid equations 
experiencing a one dimensional adiabatic compression^- 

Appendix C: Calculation of / for isothermal and 
polytropic distributions 



Let us define 



I = 



dv. 



V 



(C.I) 



For the isothermal distribution (|I02[) . the integration is 
straightforward and yields 

hso = -Pp- (C.2) 

For the polytropic distribution (|158[) , the integrals can be 
expressed in terms of Gamma functions. After simplifica- 
tion, using F(.x + I ) = xT(x), we obtain 

7 

which includes the isothermal distribution (|C2[) as a spe- 
cial case corresponding to 7 = I. 



Ipoly 



(C.3) 



Appendix D: Direct calculation of W(ix) 



Let us compute W(z) defined by Eq. (|106|> when z = ix. 
When x > 0, using Eq. (p)3"]) . we have 

r> + 00 



(B.16) W (ix) = -^=! 

V 2-7T J-t 



t 



t — ix 



^/ 2 dt. 



(D.I) 



36 As noted by van Kampen (1957), this differs from the re- 
sults of Thomson & Thomson (1933) whose considered isother- 
mal perturbations and from the results of Gross (1951) who 
considered three dimensional adiabatic perturbations. 
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Multiplying the numerator by t + ix and noting that the 
imaginary part vanishes by symmetry, we obtain 



W{ix) 



1 



i 2 



~ t2 / 2 dt. 



For any real x, we have the identity 



(D.2) 



1 f + °° t 1 



r2w 



- ±2 i 2 dt 



l - ^| -\x 



72 



1 - erf 



V2 



Therefore, for x > 0, we obtain 

When x < 0, using Eq. ([55)1. we have 
1 f +ao t 



(D.3) 



(D.4) 



W(ix) 



12-K J-oo t- ix 

which is the same as 



W(ix) 



1 



+ OO 



' / 2 dt-V2nxe x / 2 , (D.5) 



- t2 / 2 dt-^xe x2 ' 2 . (D.6) 



Using identity (|D.3[) and the fact that erf(— x) = — crf(x), 
we finally obtain for x < 0: 



(D.7) 



Comparing Eqs. (|D.4p and (|D.7|) . and noting that W(0) = 
1, we see that such expressions are valid for any real x. 

Appendix E: The case of two cold streams 

We consider a distribution function of the form 



(E.l) 



corresponding to two symmetric cold streams (T = 0) 
moving in opposite direction with velocity ±u a . We shall 
consider the stability of this distribution with respect to 
the Vlasov equation for self-gravitating systems, plasmas 
and for the HMF model. 

Self- gravitating systems: integrating by parts, the dis- 
persion relation ([52)) can be rewritten 



4ttG f f(v) 
1 + -ttt- / — a dv = 0. 



k 2 J 

For the distribution (IE. II). we obtairr 37 



1 



1 



(E.2) 



(E.3) 



The same dispersion relation is obtained if we consider 
two gaseous cold streams described by the Euler equations (cf 
Talwar & Kalra 1966). 



This is a 4th degree equation for uj with real coefficients. 
Therefore, it admits 4 complex solutions that occur in 
pairs which are complex conjugate of one another, i.e. 
^1,2 = u r ± iuji and 0*3,4 = u' r ± ic^. Accordingly, the 
system is stable for a perturbation of wavenumbcr k iff 
the four roots are real (a;, = = 0). However, it is clear 
that Eq. (|E.3[) has no solution with real u so that the dis- 
tribution (jE.l[) is always unstable (to any wavenumbcr). 
fn fact, the roots of this equation can be obtained analyt- 
ically. They are given by 



U 



2;„2 




Let us introduce the wavenumbcr 



kn = 



irGp 



1/2 



(E.4) 



(E.5) 



If k < ko, then ui 2 is real negative so that cj is purely 
imaginary. In that case, the perturbation grows without 
oscillating (ui r = 0, Wj > 0). On the other hand, if k > ko, 
then uj 2 and u are complex. In that case, the perturba- 
tion grows and oscillates (uj r ^ 0, u)i > 0). This is called 
overstability. Therefore, two cold streams are unstable for 
k < fco an d overstable for k > kg. If we consider an asym- 
metric distribution, we find that the system is always over- 
stable. 

Plasmas: in the case of plasmas, performing the trans- 
formation G — > — e 2 /m 2 , the dispersion relation be- 



come; 



:S8 



2-ire 2 p 



k' = 



m- 



1 



1 



(E.6) 
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From the graph of the function k 2 — F(ui/k), we find that 
Eq. ()E.6[) has 4 real roots if k > k c where 



Une 2 p\ 1/2 1 
\ m 2 



(E.7) 



Therefore, the distribution ([E.l[) is stable for k > k c and 
unstable for k < k c . The roots of Eq. ()E.6[) are given by 



= v 2 k' 



2ire 2 p 



1± 



Tte 2 p 



(Ei 



The solutions with the sign + correspond to u> 2 real pos- 
itive (pure oscillations). The solutions with the sign — 
correspond to to 2 real positive if k > k c (pure oscillations) 
and to uj 2 real negative if k < k c . In that case, uj is purely 
imaginary so that the perturbation grows without oscil- 
lating (w r = 0, Ui > 0). There is no overstable mode. If 
we consider an asymmetric distribution with asymmetry 
A, we find that k 2 = Aire 2 p/m 2 v 2 as for the symmetric 
case. For k < k c , the system is overstable. 



The same dispersion relation is obtained if we consider two 
gaseous cold streams described by the Euler equations (cf Kahn 
1958). 
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HMF model: for the attractive HMF model, performing 
the transformation k — > 1 and AttG — > k/2, the dispersion 
relation becomes 



kp 
T 



1 



1 



(V a + 0J) {V a - Oj) 



(E.9) 



where k is the coupling constant. The distribution (jE.lll 
is always unstable. The roots of Eq. (|E.9|) are given by 




Let us introduce the velocity 



v = [T7 



1/2 



(E.10) 



(E.ll) 



If v a < vq, then oj 2 is real negative so that uj is purely 
imaginary. In that case, the perturbation increases with- 
out oscillating (oj r = 0, Wj > 0). On the other hand, if 
v a > ^o, then uj 2 and uj arc complex. In that case, the 
perturbation increases and oscillates (uj,. ^ 0, cjj > 0). 
Therefore, two cold streams are unstable for v a < vo and 
overstable for v a > Vq. If we consider an asymmetric dis- 
tribution, we find that the system is always overstable. 

For the repulsive HMF model, performing the trans- 
formation k — > — k, the dispersion relation becomes 



1= kp 
4 



1 



1 



(v a +uj) 2 (v a -uj) 2 



(E.12) 



From the graph of the curve F(u) = 1, we find that Eq. 
(|E.12p has 4 real roots if v a > v c where 



1/2 



(E.13) 




Therefore, the distribution (jE.lj) is stable for v a > v c and 
unstable for v a < v c . The roots of Eq. (jE.12[) are given by 



(E.14) 



The solutions with the sign + correspond to uj 2 real posi- 
tive (pure oscillations) . The solutions with the sign — cor- 
respond to uj 2 real positive if v a > v c (pure oscillations) 
and to oj 2 real negative if v a < v c . In that case, uj is purely 
imaginary so that the perturbation increases without os- 
cillating (uj r — 0, uji > 0). There is no overstable mode. If 
we consider an asymmetric distribution with asymmetry 
A, we find that v 2 = (A + l)kp/A. For v < v c , the system 
is overstable. 
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